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Abstract 



We develop the hehcity modulus as a criterion for superconducting order 

in the mixed phase of a fluctuating type II superconductor. We show that 

there is a duality relation between this helicity modulus and the superfluid 

density of a system of analog 2D bosons. We show that the vortex line lattice 

exhibits a perfect Meissner effect with respect to a shearing perturbation of 

the applied magnetic field, and this becomes our criterion for "longitudinal 

superconductivity" parallel to the applied field. We present arguments based 

on the 2D boson analogy, as well as the results of numerical simulations, that 

suggest that longitudinal superconductivity can persist into the vortex line 

liquid state for systems of finite thickness, comparable to those commonly 

found in experiments. 
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I. INTRODUCTION 

The mixed state of a type II superconductor in an applied magnetic field H is charac- 
terized, in mean field theory, by a spatially varying order parameter ip{r) whose amplitude 
vanishes continuously as Tc2{H) is approached from below.H While this description is ade- 
quate for traditional superconductors, the importance of thermal fluctuations in determining 
the behavior of the high temperature superconductors^ now requires one to find a reasonable 
criterion for superconducting coherence in the mixed state, that is defined in terms of an 
average over all fluctuating configurations ip{r). One possibility is the correlation function 
{ip*{r)ip{0)). However controversy has arisen over the proper gauge invariant definition for 
this correlation function;EnZl the most straightforward definition leads in three dimensions 
(3D) to correlations which decay exponentiallyou (albeit with a long decay length) even 
within the Abrikosov vortex line lattice state, once harmonic elastic fluctuations of the vor- 
tex lines are included. The flux flow resistance of an unpinned vortex lattice in a completely 
clean material also is contrary to the conventional idea of a superconductor as a state with 
zero resistance. In this paper we propose using the helicity modulus as a clear equilibrium 
quantity that can distinguish superconducting from normal behavior in the mixed state. We 
will show that the helicity modulus, which for a neutral superfluid is proportional to the 
superfluid density,^ is in a superconductor (or charged superfluid) related to the magnetic 
susceptibility of the system to a small perturbation in applied magnetic field, about the 
uniform applied H. Recall that it is the magnetic response, rather than electrical resistivity, 
that gives the true defining equilibrium signature of the Meissner transition in either a type 
I or type II superconductor. Here we will show that for the mixed state of a type II super- 
conductor, the vortex line lattice displays (in the absence of dislocations) a perfect Meissner 
effect with respect to a certain type of shear perturbation of the applied field, for which 
the screening currents run parallel to H. Such a shear Meissner effect, also referred to as 
longitudinal superconductivity^ we will take as the defining equilibrium criterion for super- 
conducting order within the mixed state. We will then give a set of arguments, including the 



results of numerical simulations, that suggest that for system sizes of experimental interest, 
longitudinal superconductivity can persist above the vortex line lattice melting transition, 
into the vortex line liquid state. Some of our results have been briefly presented earlier ,01121 
for the case of an isotropic system. Here we provide much greater detail, and generalize our 
formalism to the uniaxial anisotropic case. 

The rest of this paper is organized as follows. In Sec. II we define our London model 
for a continuous anisotropic supercondutor, giving the mapping of the Hamiltonian from its 
representation in terms of the wavefunction phase angle, to its representation in terms of 
interacting vortex lines. In Sec. Ill we define the helicity modulus T^i,(q), discuss its relation 
to the magnetic susceptibility, and describe the important physical parameters that may be 
extracted from it. We also discuss in some detail the mapping between the interacting vortex 
lines and an analog system of interacting two dimensional (2D) bosons.llj We show that an 
interesting duality exists between the helicity modulus of the 3D superconductor model and 
the helicity modulus of the 2D analog bosons, for both the superconductor with a finite 
magnetic penetration length A (our main concern in this work), and the superconductor in 
the A ^ oo approximation. In Sec. IV we analyze the helicity modulus within the elastic 
approximationli3'y for a vortex line lattice, and demonstrate the existence of the shear 
Meissner effect. We calculate how the penetration length for the shear perturbation increases 
with temperature, due to second order elastic fluctuations. In Sec. V we consider the vortex 
line liquid, and show how the hydrodynamic approximations yields the disappearance of 
the perfect shear Meissner effect. We discuss how the Kosterlitz-Thouless transition^ of 
the analog 2D bosons can yield a cross-over from a normal vortex line liquid to a line liquid 
with longitudinal superconductivity, and estimate this cross-over temperature as a function 
of system thickness and applied magnetic field. In Sec. VI we discuss our numerical Monte 
Carlo (MC) simulations. We define the Hamiltonian and helicity modulus for a discretized 
lattice superconductorf^ini discuss our MC algorithm, and present our numerical results 
for an isotropic model. We compute the helicity modulus and other measures of vortex 
fluctuations, and find evidence for longitudinal superconductivity within the vortex line 



liquid state. In Sec. VII we present our conclusions and discussion. 

II. LONDON SUPERCONDUCTOR MODEL 

We will model our uniaxial superconductor as a three dimensional continuum with the 
weak coupling direction parallel to the z axis. The bare magnetic penetration length along 
this weak direction is A^, while A_l is the penetration length within the more strongly coupled 
xy planes, t] = A^/Aj, is the anisotropy parameter. 

The Ginzburg-Landau Hamiltonian for the Gibbs ensemble, in the London approxima- 
tion, can then be written as. 



n\e,A] = h^j d?r 



V X (A - A'^ 



(1) 



where 

Jx ^ ^^/(levr^Ai) (2) 

is the couphng within the xy plane (0o = hc/2e is the flux quantum), 

ri^=r]y = l, r]^ = r] = \JX± (3) 

deflne the anisotropy, and (0o/27r)A and (0o/27r)A'^^'* are the vector potentials for the in- 
ternal and applied magnetic flelds, 

V X A = 27rb, V X A^^* = 27rh , (4) 

where b = B/0o and h = H/0o are the densities of flux quanta of the magnetic fleld B(r) 
inside the superconductor, and the externally applied fleld H(r). In Eq.(P, 9{r) and A(r) 
are thermally fluctuating variables to be averaged over in a partition function sum, while 
^ext^j,^ is a flxed (quenched) fleld. 

It will be useful to introduce the induced magnetic vector potential, 

^ind ^ ^ _ ^cxt^ y ^ ^ind = 27r(b - h) = 27rb'°'^ , (5) 



in terms of which the Hamiltonian becomes, 



n[e,A] = ^j^ I dh 



T.%'iy,e-A^^'-j^ff + x 



Vx A 



ind 



(6) 



and the partition function is to be viewed as a sum over A™'^. 

The Hamiltonian of Eq.(^ can now be mapped onto to a model of interacting vortex 
lines. If we define the superfiuid velocity, 



V = V^ - A""* - A"^'^ = V^ - A 



(7) 



then by Eq.(HD we have 



V X V = 27r(n - h - b^'^^) = 27r(n - b) , 



(8) 



where n = ^ V x \/6 is the vortex line density, consisting of singular lines of integer vorticity 
in the phase angle 0{v). 

Defining the Fourier transforms. 



Vg = / c?V €*''■'■ vfr) 



v(r) = ^E 



-iqr 



V 



(9) 



(where V is the system volume), we can then write the Hamiltonian (^ as. 



^ = ^ E [%\,v.,, + 47r2Ai6-^6-J^; 



(10) 



9,/^ 



and solve Eq.(P) as. 



Vg = — 27ri 



qxg + 



cp^_(rv2ViK;;;'' 

q2 



(11) 



where x(r) is any smooth scalar function. Substituting Eq . ([TT|) into Eq.(0), and completing 
the square in Xqi results in, 
47r2j, 



H 



2V 



E 



q L 



^indl 



K-hg-bf]-V^"-[n_,-h_g-b':^, 



+(gi + r'qiW,^x~, + Aib--^ ■ b- 



(12) 



where 



V0 = 1 
1 q2 



I 



A^ - A_L 



[z X q)(z X q) 



(13) 



Xlql + XlqV 
is the "bare" vortex line interaction tensor, before screening by magnetic field fluctuations, 



and 5xq = Xq~ X^q is the fluctuation of Xq away from the value 



Xq 



{\i-\i)qM^{^,-\-Kn] 



(14) 



q\\W + Aig2) 

Substituting x° into Eq.(pTD gives the superfluid velocity v° that minimizes 7i for a given 
configuration of n„ — h„ — bj°*^. 5xq represents a smooth "spin- wave" fluctuation about this 



We can now complete the square in b™'^ in Eq.(|^), subject to the constraint that q-b 
0, to get, 

47rV 



ind 



n 



2V 



E 



9 L 



[rig - h,] ■ Vg • [n_q - h_g] + (g_L + r/ qz)^Xq^X-q + (^K • Ug • (5b 



(15) 



where the tensor 



^ g ^q ^q ^q ^q '^±^q ^q 



\' 



1 + Aig 



2 ^2 



I 



A?-A2 



:(z X q)(z X q) 



(16) 



1 + \lql + \lql 

is the uniaxial anisotropic generalizationEJ of the familiar London vortex line interaction, 
and 



U, = A2l + vo 



(17) 



is the interaction tensor for fluctuations of magnetic field, 5hq = b™ — b™ ' , about the 
value 



bind,o ^ u-1 . V; • K - K] = ^V, . [n, - hg] 



(18) 



that minimizes Ti. for a given configuration rig — h^. 

Eq. ( [T3| ) represents the Ginzburg-Landau Hamiltonian written in terms of decoupled spin 
wave, magnetic field, and vortex line fluctuations. The partition function is to be summed 
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over all smooth Sxq-, all smooth 5hq subject to the constraint q ■ 5hq = 0, and all singular 
vortex line distributions rig with conserved vorticity q ■ n^ = 0. 

The interaction Vg of Eq. (0) is given as a tensor, with non- vanishing off diagonal com- 
ponents. However, as shown by Carneiro et al.f^ one can exploit the conservation of vorticity 
to rewrite Vg in a purely diagonal way. Using 

ng ■ (z X q)(z X q) ■ n_g = q]_ng± ■ n_g_L - [q_L ■ ng_L][q± ■ n-q±] 

= glng_L • n_g_L - g^rig^n.g^ , (19) 

where n_|_ = {nx,ny) is the transverse part of the vorticity, and a similar result for hg, we 
can rewrite the vortex line interaction part of the Hamiltonian as, 

47r^J_L 
^v = -izTT- Y. Vgf^Kf^ - Kt^] [n-gf^ - ^-<?m] ' (20) 

where Vg^: = Vqy = Vg±, and, 

V ^i V >^l(^ + >^y) (21) 

In most of this paper we will be considering behavior in the presence of a uniform applied 
magnetic field, for which hg^o = and hence bg'^Q = bg^o- 

III. HELICITY MODULUS, MAGNETIC SUSCEPTIBILITY, AND 2D BOSONS 

A. Definition of the Helicity Modulus 

If we define the supercurrent as, 

jm = J±v;:\ = J^%\^,o - A-* - Af) , (22) 

then from Eqs.(^0|) and ([2^ ) we see that 

where J^ = —TlnZ is the total free energy for the partition function Z = j T>6T>A™'^e~^/^ . 



Consider now a small perturbation about a uniform applied magnetic field hoz, with 



A<= 



2'Khoxy + 5A^^^. We define the helicity modulus T^i,(q) as the linear response 



coefficient between the induced supercurrent and the perturbation 6 A 



ext 



{jq^.) = -T^,^(q)M; 



ext 
qu 



(24) 



From Eqs.(E3D and 



we then have 



T^M^lq) 



^{Jq^,) 



dA'^f 



V 



92^ 



9A-t9A-v 



V 



d^n 



^M-tM-V/o 



V 
T 



dn dn 



&H 



dn 



^M-t M-V/o 



ft Acxt / \ ft Acxt / 



(25) 



where the subscript "0" denotes the unperturbed system with (5A^^* = 0. For a uniform 
system, the third term on the right hand side of Eq.(|^) may be ignored as {iq^)o = 
—V{dT-C/dA'^_^g^)Q = (for the mixed state, we are assuming that q is smaller than any 
of the reciprocal lattice vectors of the vortex lattice). 

Applying Eq.(^) to the Hamiltonian ([T0|), and using the definition of Eq. (^2]) , then 
results in, 



^^/.^(q) = J±Vu 



s _ :^i!k (v V )' 



(26) 



The form of Eq. (pG]) , expressing T^,y(q) in terms of a velocity correlation, is familiar as 
defining the superfiuid density of a neutral superfiuid, or equivalently the helicity modulus 
of an XY mo del. B 

Alternatively, we could apply the results of Eqs. fpSf ) and (pSJ) to the form of TC in Eq.(|l]) 
to get. 



(j,,) = -J±Al([qx(qxA^^ 



-27r2j_LA^([q X b 



-"g i^J■/ ' 



(27) 



and. 



T^M^^lq) = -^-lAj 



q'^Sf,^ - q^qy - 



47rVxAi 
VT 



[q X bg]^[q X b_ 



q\u 



(28) 



where we have used bg = b™ for the unperturbed system, and (bg)o = for finite small 
q. Eq.(^) is just a statement of Ampere's law, relating the induced magnetic field to the 



flowing super currents. Eq.(]28|) expresses T^,y(q) in an explicitly gauge invariant form, in 
terms of correlations of the fluctuating internal magnetic field b. 

Finally, we can also express T^jy(q) in terms of vortex line correlations. Using the form 



of 7i in Eq.(|T5[), substituting in 27rhq 
as in Eg. (PSJ) , results in. 



iq X AJj , and taking the appropriate derivatives 



Tr^^(q) = Ji 



(q X z>) ■ Vg • (q X fi) 



VT 



(q X z>) • V, ■ (n_gn^) ■ V^ ■ (q x /i) 



(29) 



where Vg is the vortex line interaction tensor of either Eqs. (|T^) or (0). 

Note that the hehcity modulus is Hermitian, T^,^(q) = T*^(q) = T,^^(— q). Also note 
that any longitudinal component of ^A^^* produces no response in (jg), since T(q) ■ q = 0. 
This is as expected since a longitudinal component of (5A^^* produces no magnetic field, and 
can be eliminated by a gauge transformation. Henceforth it will be simplest to work in the 
London gauge in which q ■ A^^* = 0. 

The tensor products in Eq.(p9D can be simplified greatly if we restrict our interest to 
wavevectors lying along the symmetry directions, i.e. q = gx, gy, and gz. Changing 
notation for the sake of clarity, from n^ to n(q), and using Eq.(pTD for Vg, we find for the 
diagonal elements 

Jl^W \. 47rVxAi {n„{qi))n„{-qv))Q 



Tf,{qi)) = T^^{q9) 



1 + xw 



(30) 



VT 1 + A2g2 

where fi, z/, a are any cyclic permutation of x, y, z, and A^ is either A^ or A_l depending on 
whether ^ = z or ft = x, y. Note that T^(gz>) ~ g^ as g -^ 0. 
The off diagonal elements are 

J±XW r47r2j_LAi {n^{qa)n^{-qa))o' 



'^^^A(1^) 



1 + xw 



(31) 



VT 1 + A2g2 

However for q = qa, q- n^ = implies that n^{qa) and ny{qd) fluctuate without constraint, 
and since the free energy of Eq.(|20|) is symmetric separately in n^ -^ —n^ and in n,^ — > —riy, 
we expect that the off diagonal terms will vanish. 

Henceforth we will restrict ourselves to the simple cases given by Eq.(^). For a uniform 
applied magnetic field along the z direction, and taking here and henceforth /i, v, a to be a 



cyclic permutation of x,y,z, we have the three distinct cases, (a) Ty(gz), (6) T^^qy), and 
(c) Tz{qx.). In Fig.|l| we show a schematic of the magnetic field lines corresponding to these 
three different perturbations. As suggested by these diagrams, we will refer to (a) as the tilt 
perturbation, (6) as the compression perturbation, and (c) as the shear perturbation. We 
will find that the first two cases are determined by the transverse and longitudinal magnetic 
susceptibilities respectively. We will find that the presence of a perfect Meissner effect with 
respect to the shear perturbation is a convenient measure of superconducting coherence for 
the mixed state. Because the screening currents involved in the shear perturbation run 
parallel to the applied magnetic field, a perfect Meissner effect for the shear perturbation 
has also been termed longitudinal superconductivityB 

B. Magnetic Susceptibilities and Renormalized Penetration Lengths 

As indicated above, the helicity modulus T^,^(q) is closely related to the magnetic sus- 
ceptibility. Combining Ampere's Law (p7| ) with the definition of T^,^(q) in Eg. p^) , we 
have 

(^Af) = --^-5A-. (32) 

For the three cases of Eq.(^0|), corresponding to perturbations SAJ^^{qi)) where /i, i^, o" is 
a cyclic permutation of x,y,z, T(q) is diagonal and so we can substitute into the above 
27iSb^^'^{qiy) = -iqSA^;^'^{qu) and 2Tx5K{qi)) = -iq5 A'^^^qu) to get, 

T^(gz>) _ d{U:''{qu)) 



A-nxM^) , (33) 



J^\\q^ dK{qv) 

where h™'^ /An is the induced magnetization, and Xa{q^) is the magnetic susceptibility at 
wavevector qv for a perturbation in applied magnetic field in direction a. 

To get a feel for the information contained in the helicity modulus, or equivalently the 
magnetic susceptibility, consider first the case of zero field, in the absence of vortex line 
fluctuations. When n^ = 0, Eq.(PDD yields T^ = J^\\q^ /{I + A^g^). Substituting into 
Eq.dD gives, 

10 



4vrXa(<?'>) 



(34) 



1 + \y 

This describes a perfect Meissner effect. As q -^ 0, {5U^'^{qi')) = ATxXa{,qi^)5hc,{qu) = 
—6h(j{qu), and so the total field inside the superconductor, {6ba) = Sh^ + (56™'^) vanishes. 
The perturbation 6h„ is completely expelled from the system. The length scale on which 
this expulsion takes place is A^. 

In the presence of vortex line fluctuations, we can write a phenomenological form for the 
helicity modulus at small q, 



In this case, substituting into Eq.(|33|) one gets. 



(35) 



Airxaiqi') 



-7m 



(36) 



Now only a fraction 7^ of the applied perturbation is expelled from the system; this expulsion 
takes place on the length scale A^r. We thus see that —7^ gives the long wavelength 
magnetic susceptibility, while A^r is the magnetic penetration length, as renormalized by 
vortex fluctuations. 7^ and A^r are the important physical parameters to be extracted from 
the helicity modulus. 

Formally, we can define 7^ and A^r in terms of the small q expansion of the vortex line 
correlation that appears in Eq.(0). If we define. 



{na{qi>)n^{-qi>))o = n^o + n^ig^ + ri^2'?^ + 



(37) 



then we have, 



and. 



X2 



7^ = - limAnXaiqi^) = 1 ;^^ — n^^o , 



q^O 



lim 



dxa{qi>) 



>^l,Xa{qi>) dq^ 



VT 



47rVxAi {n^Q-n^i\^'^) 



VT 



7m 



(38) 



(39) 
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Thus 7^ = 1, or equivalently n^o = 0, signals a perfect Meissner screening of the perturbation 
6A^^^{qi)). For zero apphed magnetic field, this has a simple physical interpretation: one is 
in the Meissner state if there are no infinitely large vortex rings. 

Although the helicity modulus should have the form of Eq.(|35D both below and above the 
superconducting transition, due to the presence of ordinary fluctuation diamagnetism above 
the transition, we expect that a phase transition will be indicated by singular behavior in 
the parameters 7^ and A^r. In particular, a transition from a state with perfect a Meissner 
screening of the perturbation 6A'^^^{qi)) will be signaled by a singular decrease of 7^ from 
unity, as well, presumably, by a divergence in A^pj^. For such a case, it is reasonable to 
interpret Ug = m^c^ /Aiie^X^^^ as the density of superconducting electrons. 

We stress at this point that 7^ and A^r are describing the response of the system to a 
small spatially varying perturbation about a uniform applied field and not the response to 
this uniform field itself. 

C. 2D Boson Analogy 

Much work on vortex line fluctuations has been done utilizing an analogy between the 
magnetic field induced vortex lines in the mixed phase of a three dimensional superconductor 
and the imaginary time world lines of two dimensional bosons within a Feynman path integral 
description of 2D quantum mechanics. Here we will show the explicit connection between 
the superfluid density of these analog 2D bosons and the helicity modulus T^(gx) giving the 
response to the shear perturbation of Fig. |I]c. 

In this analogy, as introduced by Nelson,0 the energy of vortex line fluctuations is mod- 
eled by two pieces: {%) a line tension representing the vortex core energy and self interaction, 
and {ii) a pairwise interaction between all vortex line segments which lie in the same xy 
plane at equal heights z. This simplified vortex interaction is expected to be a reasonable 
approximation when the vortex lines remain, over the length scale A, approximately parallel 
to the applied field. This simplified vortex interaction is then mapped into a 2D boson 
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mass and an instantaneous pairwise boson interaction. The mapping results in the following 
correspondences (quantities on the left refer to the 2D bosons, those on the right refer to 
the 3D superconductor): imaginary time, r ^^ z, distance in direction of applied field h; 
^boson/^boson ^^ L^, length of system parallel to h; /iboson ^^ T, temperature of 3D supercon- 
ductor; boson mass, mboson <-> ci ~ vrJ^, where ei = "^"^ei, and ei is the single vortex line 
tension. 

In Appendix ^ starting from the standard definition of the superfluid density as the long 
wavelength limit of the transverse momentum susceptibility,li3c2l we derive an expression 
Eq.(^) for the number density of superfluid bosons, Ps boson, for a system of 2D interacting 
bosons, expressed in the form of a path integral over boson world lines. In Eqs. (|A9|) and 
( lAlOp we show that p^ boson is related to the helicity modulus of the 2D bosons, Tboson(Q'), 
by limg^o Tbosonl?) = (^boson/"^boson)Ps boson- Wc uow recast the resuhs of Appendix into 
the language of vortex lines. 

For a magnetic field induced vortex line i parameterized by its transverse deflection ri±{z) 
in the xy plane at height z, the vortex line density is given by, 

dri±{z) 



n{r^,z)=Y,6^'\r^-r,^{z)) 



z + 



dz 



(40) 



Using the above correspondences between the analog bosons and the superconductor, we 
then have for the term that appears in the boson path integral of Eg . i\K^) , 

/''^bosoii/ -^ boson ,^ . (i.T^^. 

I drY^'^d'^^^ = ny{q±) . (41) 

Eq.( [A10| ) for the 2D boson helicity modulus can then be written as, 

-Tf^ — = T2-K(^^)^?y(-^^))o > (42) 

-'- boson -^_L 

where L± is the length of the system in the xy plane. The vortex correlation that appears in 
Eq.(^2D above is precisely the same correlation that enters Eq.(|30|) for T z{qyi), which gives 
the response to the shear perturbation of Fig.|T]c. Taking the limit g -^ in Eq.(^) and 
combining with Eqs. (|37|) , (|38|) and ( [AlOj ) then gives, 
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2 7 \2 



7. = 1- 



477 V, A 



L-T 



Tboson(g -^ 0) 



Tv 



boson 



2 T. \2 



1- 



47rVxA 



"■bosonP- 



s boson 



^^boson-' boson 



(43) 



This leads to the following identifications, originally pointed out by Feigelman and co- 
workers:El the 2D boson superfiuid phase with Ps boson > corresponds to a 3D vortex line 
normal diamagnetic phase with 7^ < 1; the 2D boson normal fiuid phase with ps boson = 
corresponds to a 3D vortex line phase with 7^ = 1, and hence with longitudinal supercon- 
ductivity characterized by a perfect Meissner effect for shear perturbations. 

Having made the above observation, there now exists the possibility, as first suggested 
by Nelson,EJ that a Kosterlitz-Thouless (KT) transition from superfiuid to normal fiuid in 
the analog 2D boson system, could appear in the 3D superconductor as a transition from a 
normal vortex line liquid state to a vortex line liquid with longitudinal superconductivity. 
Fisher and Lee,lHHl and more recently Tauber and Nelson,B3 have argued that if one relaxes 
the periodic boundary conditions along z that is assumed in the boson analogy, and uses 
instead the free boundary conditions which are more realistic for a bulk 3D superconductor, 
the sharp KT transition no longer exists. Nevertheless, one might expect that a clear cross- 
over remnant of this KT transition should still be observable in the superconductor. We will 



return to discuss this KT cross-over in Sec. VB 



D. A ^ 00 Approximation 



Many numerical simulations,c3E2l as well as other theoretical approaches such as the 
"lowest Landau level" approximation,!^ have been based upon the approximation of taking 
Ax — »• cxD, while keeping Jj_ finite. This approximation corresponds to taking a spatially 
uniform internal magnetic field b which is equal to the applied field h. Such an approxima- 
tion can be shown to be exact for modeling the analog system of a 3D neutral (uncharged 
particles) superfiuid in a rotating bucket. Q't3 It is interesting to see how the helicity modulus 
and the 2D boson analogy look within this A — > 00 limit. 

In this case, the interaction between vortex lines is given by the "bare" interaction tensor 
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Vg of Eq.(|T3D. One can show that the correct hehcity modulus T^(gz>) is obtained by taking 
the hmit A^ ^ oo in Eg. (|30|) , keeping J^ and A_l/A^ constant, 

47rV^ {n„{qi))n„{-qv))Q 



^M^) = Jf^ 



(44) 



VT g2 

where J^ = </±(A_l/A^)2 is the couphng in direction fi. 

Noting that Ti.^ = (47r2 J_|_/2V) J2q n^ ■ V° ■ n_g must have a finite thermal average, and 
since V^ ~ l/q"^, it must therefore be true that as g ^ 0, 

{na{qi')n^{-qi>))o ~ q^ ■ (45) 

Substituting Eq.(^5D into Eq.(0) we see that, in contrast to the finite A^ case where we 
found T^(gz>) ~ g^ as g ^ 0, here we find that liniq^o ^mI?^) is in general a finite number. 
This differing dependence of the helicity modulus on q, in the small q limit, is one of the 
characteristic differences between a charged superfluid (with a finite A giving a coupling to 
a fiuctuating vector potential) and a neutral superfiuid (with A ^ oo leading to a frozen 
vector potential). For the A ^ oo case of a neutral superfiuid, T^(0) = \im.g^oTfj^{qi)) is 
just proportional to the number density of superfiuid particles, as discussed in Appendix 
^ for the two dimensional case, and used in the preceding section. We therefore expect to 
find Tz{0) > for T < Tc in an ordered phase with longitudinal superconductivity, and 
^''2(0) = for T > Tc in the normal phase. 

We now consider the 2D boson analogy for this A -^ 00 approximation. Combining 
Eg. (^2]) for the helicity modulus Tboson(Q') of the analog bosons with Eg. (^5)), we see that 

Tboson(g) ~ g^ for small q . (46) 

Thus the analog 2D bosons have a helicity modulus characteristic of a 2D charged superfiuid! 
This is in agreement with the results of Feigelman and co-workers,^ who show that the 
system of vortex lines interacting with the true London interaction of Egs.(pOD and ([2T| ) (as 
opposed to the more simplified interaction of Nelson's model) can be viewed as a system of 
analog 2D bosons whose interaction is mediated by a massive vector potential. As A — i> cxd, 
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the mass associated with this vector potential vanishes, and one has a system of 2D charged 
bosons interacting with 2D electrodynamics. 

We can develop the analogy further. Combining Eqs.(^) with ( ^^ we have, 

4:7l'^Jz Tboson(g) 



Tz{qx) = Jz 



1 - 



-'^2-' H ^ boson 

One can then define the proportionality coefficient 7boson of Eq. (^) by. 



(47) 



-I- boson I, y J /boson -^ boson . o t /bo 

47rV, 



boson 



r , (48) 



.4vr2j, 

where we have used the correspondences between superconducting variables and analog 
boson variables to arrive at the last equality. Note that when 7boson = 1 we have Tz(0) = 0, 
and when 7boson < 1 we have T ziQi) > 0. 

One can now show,EHl at least in the isotropic case, that the term [^boson/^^r^J^] = 
[^boson47rA^/0o] which appears on the right hand side of Eq.(^8D is just twice the magnetic 
energy coupling of the analog magnetic field of the 2D electrodynamics. We can rename 
this coupling [JA^jboson in analogy with the magnetic energy couphng of our original 3D 
superconductor of Eq.(|l|). Eq. (|i8|) then becomes, 

'Tboson(q') = 7boson[</A jboson? (49) 

in complete agreement with the small q limit of the form of the helicity modulus for a 
charged superfiuid, given in Eq . (^51) (as derived for our original 3D superconductor at finite 
A). — 7boson is therefore the magnetic susceptibility of the analog 2D charged bosons. To 
next order in g^ we expect, in analogy with Eq. (pSf) , that Tboson(?) has the form. 



[^A J boson Q' 

1 _|_ \2 ~2 

^ "T ^RbosonV 



ry / \ [^ ■^ Jbosony (Kf)\ 

-'-bosonl^yj 7bosonTj | 72 2 ' V / 



where Ar boson is the magnetic penetration length of the analog 2D charged bosons. If we 
take A^boson = 1/(47^^ boson), with n^boson the number density of superfiuid charged bosons, 
then combining Eqs. (|47|-pOD one can recover all the results found in Sec. VB.3 of Blatter et 
al.p^ which are therefore seen to apply in a strict sense only to the A — » oo approximation, 
rather than to the finite A case. 
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We thus have the following amusing duality. For finite A we have seen in the preced- 
ing section that the 3D superconductor, which is a charged superfluid problem, maps onto 
analog 2D bosons, which is a neutral superfiuid problem. The 3D longitudinal superconduc- 
tivity transition maps onto the 2D superfiuid transition. A perfect Meissner effect for shear 
perturbations in the 3D superconductor, with 7^ = 1, represents the normal fluid state of 
the 2D bosons with Tboson(q' ^ 0) = 0; the loss of this perfect Meissner effect, with 7^ < 1, 
corresponds to the superfluid state of the 2D bosons with Tboson(q' — > 0) > 0. For A ^ 00 
however, the 3D superconductor, which now behaves like a neutral superfluid problem, maps 
onto 2D bosons interacting with 2D electrodynamics, which is a charged superfluid problem. 
The 3D longitudinal superconductivity transition now maps onto a Meissner transition of 
a 2D superconductor. The normal state of the 3D superconductor, with T^iq — ^ 0) = 0, 
corresponds to a perfect Meissner state of the charged 2D bosons, with 7boson = 1; the 3D 
superconducting state, with Tz{q — *> 0) > 0, corresponds to the loss of this perfect Meissner 
effect for the 2D charged bosons, with 7boson < 1- 

Note that for the analog 2D charged bosons of the A ^ 00 approximation, vortices in 
the 2D condensate wavefunction will interact with a potential that decays exponentially on 
length scales greater than Ar boson, due to the screening by the 2D analog magnetic fleld. A 
vortex anti-vortex pair will therefore have a flnite energy barrier for unbinding, and so at 
any flnite Tboson there must be free vortices which will destroy the 2D Meissner state. Only 
at Tboson = {Lz — > 00), does there remain the possibility of a sharp Meissner transition 
in this 2D analog boson system, as ^boson varies. Such a transition, if it exists, is driven by 
quantum and not thermal fluctuations and so it is not in the Kosterlitz-Thouless universality 
class. We believe that it is this transition at Tboson = 0, in the A -^ 00 model, that the 
work of Feigelman and co-workers" pertains to. Recently, Tesanovic has arguedl^ that such 
a transition must be driven by the proliferation of closed vortex rings (boson anti-boson 
virtual pairs), which are left out of the naive 2D boson mapping, and that the transition 
will be in the universality class of the ordinary 3D XY model. 

The above considerations suggest that taking the A — > 00 limit in our model is rather 

17 



subtle and possibly leads to discontinuous changes in the phase diagram, although any such 
discontinuities will likely be obscured in a finite size system by very strong cross-over effects. 

IV. THE VORTEX LINE LATTICE: ELASTIC APPROXIMATION 

We consider now the mixed state of a type II superconductor. At low temperatures, such 
a state is described by the familiar Abrikosov vortex line lattice. In this case, we can evaluate 
the vortex line correlations that appear in the expression for the helicity modulus by using 
the well known elastic approximation.ll3'E3 It is now convenient to work in the Helmholtz 
ensemble, with a fixed uniform density bo of magnetic field induced vortex lines, 

6„ = f. a.= (|)"'^ (51) 

where a^ is the lattice spacing between lines in their ground state triangular lattice. We will 
denote thermal averages in this ensemble by (...), dropping the subscript "0" that we used 
earlier. 

In the elastic approximation, one assumes that vortex line excitations consist only of 
fiuctuations of the magnetic field induced vortex lines, transverse to the direction of the 
uniform apphed field. Such fiuctuations are described by the displacement field Uj(z), which 
gives the transverse displacement in the xy plane at height z, of the vortex line away from 
its position Rj in the ground state vortex lattice. The vortex line density is thus given by 
Eq.(^OD, making the substitution rj_|_(z) = Rj + Uj(z). 

If we define the Fourier transforms 

U" = r f dzY^e''^'^^^+^--''^'>u,{z), u,{z) = l^e-^('?^^+''-^')u, , (52) 

bo J i V q 



'■q 



where q^ = {<ix-,<iy) and the sum over q^ is restricted to the first Brillouin Zone of the 
Abrikosov lattice, then to lowest order in u the vortex line density at small finite q may be 
written as 

rig = iho[q, ■ UgZ - q^n^] . (53) 



Substituting the expansion for n^ in terms of the u^ into Eq.(^), summing over reciprocal 
lattice vectors, and keeping only terms up to order u^, results in the free energy functional 
for elastic vortex line displacements, 



ijo/3 



2V 



Yl \ [c44(q)g^ + Cu{(i)ql]ugLU-gL + [C44(q)g^ + C6e{(l)ql]ugTU-gT 



(54) 



Here UgL = q ■ Ug is the longitudinal part of the displacement, and UgT = \ug — qUgil is 
the transverse part. C44(q), Cgelq), and Cii(q) are the wavevector dependent tilt, shear, and 
compression elastic moduli respectively. We can now use Eqs. (^3|) and (|^) to evaluate the 
vortex correlations of the helicity modulus of Eg. (pOD , for the three simple cases of Fig. |I[ 



A. Tilt Perturbation 

We first consider Ty{qz) which gives the response to the tilt perturbation of Fig. |l|a. 
Using Eq.(^), the relevant vortex correlation, to lowest order in Ug, is 

b'^VT 

{n^{qz)n^{-qz)) = q%l{uT{qz)uT{-qi)) = ° .. , (55) 

C44[qZ) 

where we have used T^ei of Eq.(p4D to evaluate the displacement correlation. Expanding in 
g^ we have, 

(56) 



n. 



blVT 



Hyl = 



hlVT dcM 



clM dql 



<?=o 



'' C44(0) 

Combining Eqs.(|^), (|38D , (|39|) , and (^) with ( ^61) above then gives for the helicity 
modulus parameters, 

52 



ly 



47rc44(0) 



(57) 



and. 



\2 



5^ 



(iC44 



47rc44(0)7j^ 



C44(0)A; 



(58) 



<7=0_ 
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Note that from general thermodynamic argumenta23 one has, 

C44(0) 



A'K dB\ 



(59) 



where the dH±/dB± is evaluated at the average magnetic field Bqz. Hence 7^ is determined 
by the transverse magnetic susceptibility, 

dBi 



^' ^ dH^ ' 

as expected from our discussion in Sec. |1IIB| . 

Using our explicit results for C44 from Appendix M, we have. 



(60) 



87rA^5 



f] In 



Hc2 



47rA^fi 



^ 1 for large B, 



A I 3> a^ 



1 



'^y-\ 1 7^^ '?l^ ,TT ^/T^\ rrrr ~ ^ ^^"^ intermediate B, A^ < Ov < A, 



"V ^^ ^'z 



SttA^E 



r/-2 [2 ln(r7K) - 1] + 1 



1 for small B, 



Xz <^ Ov 



(61) 



For Aj^R, using Eq.([55|) and our results for C44 from Appendix ^, we have for large 
magnetic fields, A^ ^ Ov, 



r 1 



X2 



2AiA:2 



for strong anisotropy, ,„ , „ ^ r) In 

2A^H 



-2 i„ / ^c2 
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(62) 



2AiA;2 



r/ ^ In 



i/, 



c2 



for weak anisotropy 



1 



-2 i„ / ^c2 



< ?7 ^In 



5 7 ""' 2AiA;^ ' V 5 

For intermediate magnetic fields, A_l <^ Ov <^ A^, (where strong anisotropy is by definition 
implied) we have 



A 



j/R 



Ai 



=. 1 - A^ B 



'±"^0 ' 



(63) 



and for weak magnetic fields, A^ ^ Ov, we have 
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1 — 2\\kl for strong anisotropy, | ^ ^ ^InrjK 

(64) 
"' for weak anisotropy, | ^ t]^^ In r]K 



where /cg = 47ri?/0o ~ l/«v- 

For strong magnetic fields, Aj_ ^ a^, Eg . (|6^) gives Aj,r ~ l/-\/2A;o ~ ctv, independent of 
the bare Aj^. Since our definition of Xy^ in Eg . (|35D was based on an expansion in small q, it 
is doubtful that we should take such a small Aj,r too seriously as a screening length, without 
considering higher terms in an expansion in q, as well as considering the response jq+x to 
the perturbation A^^* (where K is a reciprocal lattice vector of the vortex lattice). 

For weak magnetic fields, A^ -C a^, Egs. (|63|) and (|64D give, A^r ~ A^. 

B. Compression Perturbation 

We next consider Txiqy) which gives the response to the compression perturbation of 
Fig.|I|j. Using Eg . (|53D , the relevant vortex correlation, to lowest order in Ug, is 

b'^VT 

{n,{qy)n,{-qy)) = q\{uL{qy)uL{-qy)) = °. .. , (65) 

where we have used Eg.(0) to evaluate the displacement correlation. 
We therefore have, 

■"^'-^)- '''*" 

The compression modulus in the vortex line lattice can be written as cn = cl + cqq, where 
cl is the "bulk modulus" for an isotropic compression. General thermodynamic arguments 
give,E3 

. . B dH-, , _. 



Noting from Egs. ( [Bib] ) and ( |B23| ) that, for large A_l ^ Ov, C66(0) ^ Cn, we have Cii(O) 
cl(0), and so, 
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Ix 



dH, 



(68) 



is determined by the longitudinal magnetic susceptibility. For A^ <^ a^, Eqs.( |B2^ ) and 
(|B2Bp give cii = |cl, and so 7x = 1 - l{dBjdH,). 

Using our explicit results for cn from Appendix ^, we find, 

00 



Ix 



IQtt\\B 
IGv^A^E /A 



for large B, \\^ a 



(69) 



/A \^/^ 
— e'^-/^^ for small B, A^ < 



7a, < implies that the magnetic field induced in the material is larger than the applied 
perturbation, so there is negative screening. This may be understood from Eq. ( |55| ) by noting 
that in the mixed phase one always has dBz/dHz > 1. 
For the screening length, we find 



X2 



4Xlki 



5al 
72A2 



for large B, A^ ^ a^ 



for small B, \± <^ a^ 



(70) 



Since k^ = A-nB /(pQ ~ l/ctv) both cases give X^ji ~ ia^- It is tempting to interpret this 
imaginary A^jr as indicating the rearrangement of vortex lines on the length scale a^ due 
to the penetration of the applied field, with no "healing" length at all at the surface of the 
sample. However our cautionary remarks following Eq. (|5^ , concerning the applicability of 
our results on the length scale Ov, should again be noted. 



C. Shear Perturbation 

The preceding two cases of the tilt and the compression perturbations gave information 
about the transverse and longitudinal magnetic susceptibilities, via 7^, and 7^;. However, 
since for large B A^r, A^r ~ a^ is independent of the bare A_l, it is unclear whether they 
give any information about the density of superconducting electrons, or whether they can 
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be expected to diverge at the superconducting to normal transition. A more interesting 
case is therefore given by the third possibihty, T2(gx), which gives the response to the shear 
perturbation of Fig. |I]c. 

The relevant vortex correlation we need to compute is {ny{qx.)ny{—qx.)) , but to lowest 
order in u^, Eq.(|5BD shows this to vanish identically. This is merely an artifact of our Fourier 
transform of the displacement field u, which prohibits vortex lines from having a net tilt 
away from the z axis. To avoid this difficulty we can evaluate the correlation at q = g^^x+g^z, 
with finite qz, and then take the limit as qz -^ 0. From Eqs.(^3|) and (|^) we get 

lim(nj/(q)n„(-q)) = lim — , , ^„ ° -— ^ . (71) 

g^^o' y^^^ ^" 9.-0c66(q)g2 + c44(q)g2 

For the case of a vortex line lattice with cge > 0, taking g^ — * results in a vanishing of the 
vortex correlation. From Eqs. (pH]) and ( p9D we then have A^r = A^, and 7^ = 1. We will 
see in the next section that higher order elastic corrections lead to a temperature dependent 
increase in A^r, but do not change 7^ = 1. The vortex line lattice thus exhibits longitudinal 
superconductivity, with a perfect Meissner effect for shear perturbations. 

If one assumes Cge = 0, as might describe the case of a vortex line liquid, Eq.([7TD 
results in {ny{qk)ny{—qk)) = b^VT / c^i^qyi) . In this case we have 7^ = 1 — B"^ / Attc^^^O) = 
1 — dBj_/dH± < 1, exactly as in the case of the tilt perturbation, Eqs.(|^) and (|60|). To 
summarize, 

{ny{qy.)ny{-qk)) = , 72 = 1 if Cee > 

I ( ^^ ( M\ ^0^^ 1 dB^ _ (72) 

{ny{q^)ny[-q:ii.)) = — — - , 7^ = 1 - ttt < ^ if Cge = 

C44(^gxj ari± 

The above arguments suggest that a singular decrease of 7^ from unity (or equivalently 

the singular increase of ra^o from zero), marking the loss of longitudinal superconductivity, 

serves as a convenient criterion for the superconducting to normal transition in the mixed 

state of a type II superconductor. This is one of the main results of our paper. If this 

transition is second order, we expect that A^r will diverge at the transition with Ug ~ 1/A^j^ 

the density of superconducting electrons. 
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In considering the vortex correlation at q = g^x + g^z, it is the relative order in which 
one takes q^ and q^ to zero that distinguishes between the shear and the tilt perturbation. It 
is the order corresponding to the shear perturbation that is related to the superfiuid density 
of the 2D analog bosons, Ps boson- Note that if we had defined Ps boson in terms of the q = 
winding number of Pollock and Ceperley,E3 rather than in terms of the q ^ transverse 
momentum correlation function as we have done in Appendix 0, the connection of Ps boson 
to the shear, as opposed to the tilt, perturbation would become ambiguous. 

The preceding discussion has been based upon elastic fluctuations about a perfect dislo- 
cation free vortex line lattice. Recently, Frey et alB^ have argued that, at sufficiently high 
magnetic field in a layered superconductor, the proliferation of dislocations can result in the 
loss of longitudinal superconductivity even in the vortex line lattice state. 

D. Second Order Corrections 

Our analysis in the preceding sections is based on Eq.(^), which gives an expansion of 
the vortex line density n to linear order in the displacement field u. In this section, we 
consider the effect of higher orders, by continuing the expansion in u, 

rig = ibo[q ■ UgZ - q^Uq] 



Vt7 

bo 



ib( 



2^^ .W' 



Yl {Ih- ^9'] [q ■ V'] [q • Uq-n'-n"\Z 



'■q-q'-q" 



{Qz -q'z- g") [q ■ u,/] [q • u^.] Ug_g/_g// } + . . . (73) 



for small but finite q. 

To systematically evaluate vortex correlations using this higher order expansion, one 
should also in principle extend the elastic energy of Eq. (Q) to higher order in u by taking 
the expansion above, and substituting into the vortex-vortex interaction Hamiltonian of 
Eq. (pOD . The resulting expression is rather complex. For simplicity, we will instead continue 
to use the quadratic elastic energy of Eq.(B5), however we now view the elastic moduli as 
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appropriately redefined temperature dependent parameters, in the spirit of a self-consistent 
phonon approximation. 

We consider here only the case of T2(gx), corresponding to the shear perturbation, for 
which we need the correlation, \imq^^o{ny{q)ny{—q)), with q = gx+g^z. As we have already 
seen in the preceding section, for Cge > 0, the contribution to O^u"^) vanishes. By symmetry, 
the next leading term is 0{u'^). Using the expansion of Eg . (|73D , and factorizing the average 
of the product of the four m's into all possible pairs, we find. 



,(q)n,(-q)) = blT'q^q, ^ {(g. - g;)^$-^(q')$-,^(q " q') 

+ (Q. - Q'MKyi^)^-pyi^ - q') - Ql%yi^)K'pi^)} , (74) 



where summation over a,/3 = x,y, is implied, and $ is the elasticity tensor of Eg. (^) . 
Taking Qz — > 0, keeping only terms of 0{q^), and using the fact that $ is symmetric in q as 
well as its indices, we get. 



(n,(gx)n,(-gx)) = hlT'Vq'l , (75) 



where I is the integral. 



1 i-2 1 jL2 

I = -Y " = - V - (76) 

The correlation of Eg.([75D vanishes as q^ for g ^ 0. We therefore continue to find, as in the 
preceding section, that n^o = and 7^ = 1 giving a perfect Meissner screening of the shear 
perturbation. It is straightforward to see that this result persists to all orders in u. 

However, in contrast to the preceding section, we now find a finite renomalization of the 
penetration length. Comparing the expansion of Eg.(|37D with the result of Eg.(|7|), we get 
Tt-zi = b^T'^VI. Using this in Eg. (p^) then gives. 

Thus the olu"^) term generates an 0{T) correction to A^r. Continuing the elastic expansion 
in powers of u will generate corrections to A^r in the form of a power series in T. 
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To estimate the magnitude of the correction to A^r of Eq.(|77[) we can evaluate the integral 
I using a crude approximation. The elastic moduli which appear in / are functions of 
wavevector k. However the dominant contributions to the integral will come at wavevectors 
kz ^ Y (C66/C44) k±_ and kz — Jicu/ca) k±, both giving k^ — r]k±. The dominant k± will be 
A;_i_ ~ fco = \/4vri?/0o, at the edge of the Brillouin zone. We will therefore approximate the 
elastic moduli by their values at this dominant wavevector, denoting these values as C44, Cn, 
and Cee- Within this approximation one can explicitly calculate the integral to get, 

47r C44 VC44 [Vcu + V C66 j 
Within the same crude approximation we can estimate the vortex lattice melting tem- 
perature Tm. Using the Lindemann criterionoQ'Ea that melting occurs when (m^) ~ cj^a^ 
(cl ~ 0.15 is the Lindemann parameter), and keeping only the transverse fluctuations as 
the dominant soft mode, we get, 

(79) 
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± rr 



47rCL<Vc66C44 
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Combining Eqs. ([77|), ([78|), and (|7|), and using the estimate that cn ~ cee at the 
Brillouin zone boundary, we find, 
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Using the result of FisherEa for C44 at the zone boundary, for large magnetic fields A_l ^ Ov, 
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we get. 
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(81) 



(82) 



Taking for example B = 0.2Hc2, and rj'"^ <^ 1, we estimate a 12% increase in A^j^ over A^ 
near melting, due to lowest order elastic fluctuations. Since the elastic moduli Caa(M) are in 
general larger than the values Caa, this is an overestimate. The above estimate does not of 
course include the effects of critical fluctuations near a phase transition, which for a second 
order transition should result in a divergence of A^r at Tc. 
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V. THE VORTEX LINE LIQUID 

In the preceding section we considered the vortex hne lattice at T < T^. In particular we 
showed how the response to the shear perturbation, given by Tz{qSi.), gives a useful criterion 
for superconducting phase coherence: jz = linig^o[T"z(Q'x)/J^A5^g^] = 1, or equivalently 
rizo = lim.g^o{ny{qx.)ny{—qx.)) = 0, indicates the presence of longitudinal superconductivity. 
In this section we consider behavior in the vortex line liquid at T > Tm. As a measure of 
superconductivity we focus on the behavior of Uzq. 

A. Hydrodynamic Approximation 

The simplest approximation one can make at high T is to take the Hamiltonian (^0[) and 
regard the Fourier components of the vortex line density rig as continuous, independently 
fluctuating variables, subject to the constraint q ■ rig = 0. We refer to this as a hydrody- 
namic approximation. El This yields {ny{qx.)ny{—q5i.)) = TV/[47r^J_|_Vj_(gx)]. Substituting 
into Eq.(^), we find that T^(gx) = for all values of q, reflecting the fact that well above 
the superconducting transition, an applied magnetic field will induce no supercurrents at 
all. 

For temperatures closer to, but still above T^, we expect the system to show a finite 
fluctuation diamagnetism. A better approximation can be obtained by coarse graining the 
Hamiltonian (^) over a length scale of order the inter- vortex spacing a^, and then applying 
the hydrodynamic approximation to average over the resulting coarse grained vortex density. 
This coarse grained free energy has been given by MarchettitJ as, 

^[n] = ^^ Yl {cL{<l)SnqJn-qz + C44(q)ngj. ■ n_g^} , (83) 

where Sn^ = riz — bQ, cl is the bulk modulus, and C44 is a tilt modulus of the same form as 
for the vortex lattice. Using this form we find, 

b'^VT 
Uzo = hm{ny{q±)ny{-qx)) = ^— - > , (84) 
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and from Eqs. (pH]) and (|5 



B dBj_ , . 

Thus, within this hydrodynamic approximation, the longitudinal superconductivity found in 
the vortex line lattice is lost for the vortex line liquid. Note that since Eq . (^Sf) gives 7^ strictly 
less than unity in the vortex line liquid, while 7^ = 1 in the vortex line lattice, 7^ presumably 
takes a discontinuous jump at the transition where longitudinal superconductivity is lost. 

A more detailed calculation of vortex correlations, averaging over unbounded dislocation 
loops within a continuum elastic model, has been carried out by Marchetti and Nelsoro as 
a model for a hexatic vortex line liquid. In the limit q -^ 0, the result of Eq.(p4D is again 
obtained. 

Eqs.(plD and (p5| ) are identical to the result we found in Eq.(|72D by simply taking cee = 
in the elastic approximation for the vortex line lattice. It is interesting to speculate about 
the behavior of a "soft" vortex line lattice in which the long wavelength shear modulus 
vanishes, C66(q = 0) = 0, but in which a finite shear stiffness remains on shorter length 
scales, C66(q±, g^ = 0) > for qj_ > 0. In this case, taking the limit q ^ as in Eq.(^), we 
find that Uzo = and longitudinal superconductivity remains.Ei As Marchetti and Nelsoncll 
show however, it is not possible to describe an entangled vortex line liquid with such a simple 
elastic description.^ 

B. Kosterliz-Thouless Transition 



In Sec. prn we discussed how the KT superfluid transition of the analog 2D bosons 



could appear in the 3D superconductor as a strong cross-over to a vortex line liquid state 
with longitudinal superconductivity. In his original model, NelsonEJ interpreted this KT 
transition in terms of a transition from an "entangled" to a "disentangled" vortex line 
liquid, for sufficiently thin samples. In this section, we estimate the temperature Tc for this 
KT transition as a function of sample thickness Lz and magnetic field B, and compare this 
estimate with Nelson's entanglement criterion. 
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The 2D KT superfluid transition is characterized^ by the fact that, exactly at the 
transition, the boson hehcity modulus Tboson('? -^ 0) takes a discontinuous jump to zero 
from the universal finite value Tboson/Tboson = 2/7r. Tboson is given by the vortex correlation 
of Eg. p2|) , which for a vortex line liquid can be related to the tilt modulus C44 by Eg. (p^) . 
Using Eq.(^) for C44 and applying the universal jump criterion then gives for the KT 
transition, 

^ 2 C44W ^ J^dH^ ...^ 

Thus as the thickness Lz increases, Tc decreases. In order to observe a vortex line liquid 
with longitudinal superconductivity we need the system to be thin enough that Tc > Tm. If 
we define the length, 

then we can rewrite Eq.(p6D as, 

(88) 



Tc dH± A -^2 max 



Tin dB± Lz Lz 

We thus will have Tc > Tm provided, 

dH^ , /47rc44(0)\ , 
Lz < Lz^^ = ^A = [-^-^j A • (89) 

Assuming that 044(0) in the line liquid is not too different from 044(0) in the line lattice, 
we can use our results from Appendix |B| to evaluate the length Lzmax- For large applied 
magnetic fields, such that a^ <ti \± or equivalently Hd ^ B, we have dH^/dB^ ~ 1 and so 
to leading order, 

^.max = A . (90) 

For small magnetic fields, such that Xz ^ a^, we have, 

Lzm.. = ^^^ {v'' [2 HV^) - 1] + 1} A. (91) 
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For an anisotropic material in intermediate magnetic fields, such that X± <^ a^ '^ Xz, we 
have to leading order, 

L.^.. = ^-^ {v-' [HHc2/B) - 1] + 1} A . (92) 



XIB 

Note that for a melting temperature of T^ ~ 90°K, as in YBCO, one has for large B, 
Lzmax = A ~ 1400/im. This is much thicker than typical experimental samples, which are 
generally of the order 50/im. As B decreases, L^max only gets larger. 

The above results may be compared with the original criterion for 2D boson superfluidity 
given by NelsonEJ in terms of the "entanglement length," 

e. = ^i^ • (93) 

^ 2TB ^ ' 

C,z is the length required for a vortex line to have a transverse deflection equal to the average 
spacing between vortex lines, a^. Only when C,z < Lz can the vortex lines have enough 
transverse wandering so that they may become geometrically entangled. The cross-over T^ 
between a disentangled and an entangled vortex line liquid is thus given by, 

"" 2BLz 2'K^Lz IQXIB ' ^ ' 

where we have used ei = eir/"^ with ei = ((/)Q/47rA_L)^ In k for small B. One will have Tx > T^ 
only for Lz < L'^^^^, where. 



L' =1 

z max „ 

TT 



TT^ei 



A = ^2!L!M!!)a. (95) 



^oB_ 

Except for some numerical factors, L'^^^y, agrees with L^max of Eqs. (|91|) and (0), decreasing 
as 1/B, for increasing magnetic field. In the large field limit however, our result in Eq.(pOD 
saturates to the finite value A instead of continuing to decrease. The difference between the 
results of Eqs. (pOD and ( p5| ) arises because the later is based on the the line tension ii for 



the energy of a single vortex line tilting, while the former is based on the tilt modulus 044(0) 
for the collective tilting of all lines. This points out an important distinction: geometric 
vortex line entanglement, i.e. the local wrapping of lines around each other, is not necessarily 
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equivalent to the global vortex line winding that characterizes the analog 2D boson superfluid 
phase.Ej Nelson's entanglement length of Eq.(p5D nevertheless remains the important length 
scale for local geometric entanglement, which still can have a significant effect on the dynamic 
behavior of the vortex line liquid if the barriers for vortex line cutting are high.^! 

The above discussion has been based on the familiar KT transition of an ordinary 2D 
superfluid, and predicts that as L^ -^ oo at fixed T > T^ (Tboson ^ at fixed /iboson) Tc ^ 
and so one is always in the boson superfluid state, corresponding to a normal vortex line 
liquid. FeigelmanE^ and co-workerscll however have argued that for A — > oo, the long range 
nature of the effective 2D boson interaction can lead to a normal boson fluid, and hence to 
a vortex line liquid with longitudinal superconductivity, even in the L^ ^ oo (Tboson -^ 0) 
limit for T < T^, where T^ (i.e. the critical ^boson in the boson variables) gives the 2D 
Meissner transition of the analog 2D charged bosons of the A — * oo approximation, as 
discussed in Sec. [lllEl| . Such a A — i> oo transition would probably lead to strong cross-over 



effects in the finite A case, which would obscure the KT transition when L^ ^ -^zmax; where 
the Tc of Eq. ( ]5B| ) can be very much larger than T^ . 

Searching for longitudinal superconductivity within the vortex line liquid will be one of 
the main objectives of our numerical investigations, to be discussed in the following section. 

VI. NUMERICAL SIMULATIONS 

In this section we report the results of numerical Monte Carlo simulations we have carried 
out in order to study the behavior of the system of fluctuating vortex lines. 

A. Lattice Superconductor 

To carry out numerical simulations, our first step will be to discretize the continuum to 
a cubic grid of A^ = N]_Nz sites i. The grid spacing in direction fi is taken to be 
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a^ = { (96) 

a^ = d, fi = z 

The grid spacing C,i_ in the xy plane is meant to approximate the core radius of a vortex, 
while the spacing d along z is meant to simulate the spacing between CuO planes of a layered 
high Tc superconductor. If one wants to model an anisotropic continuum superconductor, 
such as in Sec. y, one should take d = ^z, where the anisotropic Ginzburg Landau free energy 
functional givesS ^^ = i]''^^j_. Discretization of Eq.(|I|) then leads to the lattice superconductor 
model,0§ 

n[e,, Am] = E [Wm(^^+A - 0^ - Am) + '27r^C^{h^ - K^y] , (97) 

where 6i is the phase angle on grid site i, 

A,= r'A-di (98) 

Ji 

is the integral of the total magnetic vector potential across the bond at site i in direction fi, 
and if /i, z/, a is a cyclic permutation of x, y, z, then, 

2nbi^ = Ai+o^„ - Ai„ - Ai+^^^ + Aiy (99) 

is the sum of the Aj^, going counterclockwise around the plaquette at site i in direction 
/i, and gives 27r times the flux of total magnetic field through the plaquette (see Fig.||); a 
similar relation defines /ij^ in terms of Af^ . The kinetic energy piece is expressed in terms 
of the Villain function,c3 



m=— oo 



with couplings. 



a^l 




a^a„ J± = J^d, f^ = x,y 

^M = J±^^ = { ^2 (101) 



H = z 



The couplings of the magnetic energy piece are, 

32 



^n — J±^± 



c, 



X 



A 



A 



c. ^ J. - 



fj, = z 



(102) 



To express the Hamiltonian in terms of vortex line variables, we first perform a standard 
duality transformation^] of the kinetic energy piece, and then, following Carneiroll3 (in com- 



plete analogy with Egs-dlSl) — (|2ll) ), complete the square in 6; 



ind 



hq^ subject to the 



constraint that b™ is divergenceless. Our lattice Fourier transforms are defined by. 



Jqfi 



E' 



,«qri 



''^^l 1 



-'IfJ, 



^Ee— 'ft 



<?M ' 



(103) 



and the constraint that 6™ is divergenceless can be written as, Q* ■ b™ = 0, where 

g^ = 1 - e'^-^- . 
The vortex part of the resulting Hamiltonian is. 



(104) 



?^v 



2N 



/ , ' qay^qa i^qaly^—qa '^—qa\ ; 



(105) 



q,a 



where Vg^ = Vqy = Vq±, and. 



K 



9-L 



m" 



i + (^Y\Qz? + (kY\Q^ 



?J 



(106) 



K 



i + i^fm^ + ityiQA 



qz 



i + (^)'ig.p + (t^)'iQ±p 



1 + 



\Qz? + [t) IQj 



(107) 



n. 



<?M 



is the Fourier transform of the vorticity n^^ piercing plaquette i/i. Eqs. ( |105H107|) are the 



lattice equivalents of the continuum Eqs.(^) and (|2T|). 

Note that Ti^/T depends on four dimensionless parameters, which may be taken to be 
Jl/T, r] = Xz/^±, /t = A^/^j_, and ( = C,±/d. The ratio the of couplings that appear in 
the Villain kinetic energy terms of Eq. (|100| ) is then Jzr]~'^/J± = (Ax^x/A^rf)^ = {ChY- If 
one wants to model an anisotropic continuum, with d = C,z = V~^^±y then one has ( = r) 
and there are only three dimensionless parameters, J±/T, rj, and n, with JzV~'^/J± = 1- 
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Both cases are in general different from an earlier derivation of the London lattice vortex 
line interactionO which assumed equal grid spacings in all directions, a^ = ao for all /x, and 
so with C = 1 involves only the three dimensionless parameters, J_l/T, rj, and X±/aQ, but 
with JzV~'^/Ji- = V~^- Keeping the distinction az 7^ a± (i.e. d 7^ ,^_l) is essential to correctly 
model the effects of the anisotropic vortex core energy in either a continuum or a layered 
anisotropic superconductor. 

We can now define the helicity modulus for the lattice superconductor in complete cor- 
respondence with the continuum Eqs. (p2| -p^). The only change needed is to replace the 
system volume V with the number of grid sites A^, due to the slightly differing definitions 
of the Fourier transform in the continuum, Eq.(P), and on the lattice, Eq.( |103| ). As in the 
continuum we restrict ourselves to the three special perturbations of Fig.[l|, Aj^^lqi)), where 
fi, z/, a are a cyclic permutation of x, y, z. Taking the Fourier transform of Eq. (pQ]) , we get 
2'Khfj{qi>) = QlA^^^lqu). Substituting for h^ in terms of A"^^ in H^ of Eq. (|105|) , and then 
applying the definition of helicity modulus in Eq.(p5D, we get for the diagonal part of the 
helicity modulus tensor. 



^f^iQi") 



■J.[t) \Q 



and for the off diagonal part. 






TN 



(108) 



^iiu{.q^) 



l+fe) \Qr^ 



^"^ J^ [t) {n.{qa)n,{-qa))o 



NT 



1+ t \Q^ 



(109) 



Eqs. (|108| ) and ( |101J| ) are the lattice equivalents of the continuum Eqs.(^D|) and (^). The 
primary difference between continuum and lattice expressions is the substitution, g^ — >■ 
IQm^mP = (2-2cosg^a^)/a^. 

Expanding the vortex correlation at small q, 



{n„{qi))n„{-qh'))o = n^o + n^i\Q^\^ + nf,2\Qu\^ + ■■■ 



(110) 



we can again write the diagonal part as. 
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12 
'u\ 



^^^{(l^)=l^. 77^72 ' (m) 



|2 



where, analogous to Eqs. (pS]) and 



7. = 1 ^f^-.o (112) 



and 



2 /1^2 T /A 



A.rV_i 4vrV,(t^-^ 



^A«o - n^i ■ ' 



;ii3) 



Noting that Q,, ~ —iq^a^, for small g, that A^ = V/(^^(i), and that there is a slight 
distinction between Fourier components defined on the lattice versus in the continuum, 
^lattice _ ^-ij^ continuum^ g^jj ^]-^g above expressious agree completely with their continuum 
counterparts, in the limit g — > 0. 

B. Monte Carlo Method and Parameters 

To carry out Monte Carlo simulations of the lattice superconductor model, we start with 
a fixed density bo = ^^-B/0o of magnetic field induced straight vortex lines, parallel to the 
z axis. Following Carneiro, Cavalcanti, and Garter,Ej we update the system, heating from 
the ground state, by adding elementary closed vortex rings that surround only a single bond 
of the discretizing grid (i.e. rings of area ^^ in the xy plane, or area C,±d in the xz or yz 
planes). The rings are added one at a time, at random positions with random orientations, 
and then accepted or rejected according to the standard Metropolis algorithm. When a 
side of such a ring coincides with, and is oppositely oriented to, a segment of one of the 
initial magnetic field induced vortex lines, these two segments will cancel out resulting in 
a net fiuctuation of the vortex line. This procedure provides for a complete sampling of 
phase space for the vortex variables rij^ which are subject to the constraints that vorticity is 
locally conserved, J2fi[^ifi — "^j-A.^*] ~ 0' ^^^ ^^^^ ^^^ average internal magnetic field is kept 
constant, (1/A^) J2i^i = ^o^ (i.e. we are using the Helmholtz ensemble). 
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Our simulation uses periodic boundary conditions in all directions. The periodicity along 
z makes our simulation map exactly onto the 2D boson problem. In order to compute energy 
changes for the Metropolis acceptance test, it is convenient to use,E3 

AE = 27r2J^^Fi^An,^ , (114) 

where An^^ is the change in vorticity due to the vortex ring excitation, and 

Fi^ = Y, V^Xi-j - ri)nj^ (115) 

j 

represents the "potential" field of all other vortices. V^(rj) = {l/N)J2q^^^'^'^'Vqii is the 

Fourier transform of the vortex line interaction of Eqs. (|106|) and ( |107| ), where the sum is 

over all q satisfying periodic boundary conditions, g^ = 27r£^/A^^a^, £^ = 0, 1, . . . , A^^ — 1. 

In this way, the computation of AE is a local computation, involving only the sites of the 

elementary vortex ring excitation. Only when an excitation is accepted is it necessary to 

update the potentials Fj^, a calculation of order A^. Since acceptance rates are generally low 

below the transition, this method is significantly faster than a direct computation involving 

the long range vortex interactions. 

For simplicity, we have only simulated the completely isotropic case with A = A^ = A^, 
aQ = ^± = d, and hence J^ = J± = J^. Hence forth all lengths will be measured in units of 
the grid spacing oq = 1, and temperatures in units of Jq. Our simulations are for the fixed 
vortex density h^ = 1/15 whose ground state, shown in Fig.^, is a close approximation to a 
perfect triangular lattice with sides of length vTS x vT7 x vTf. We choose k = X/uq = 5, 
comparable to the vortex line spacing a^/ao ^ 1/v^ = 2.87. We study system sizes A^^ = 30 
in the xy plane, and A^^ = 15 and 30 parallel to the applied magnetic field. 

Our simulations are carried out heating from the ground state. At each temperature we 
use typically 5000 sweeps to equilibrate, followed by another 8 — 16000 sweeps to compute 
averages. Each "sweep" refers to A^ = N]_Nz attempts to add an elementary vortex ring. 
Statistical errors are estimated using the standard block averaging method. 
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C. Results: Helicity Modulus 



In an earlier reporla we presented an analysis of our data based on Eqs. ( 110 - 113 ) 



fitting our computed correlations {n^{qi))n^{—qi)))Q to an expansion in IQi/P- Here we take 
a different approach. Plotting JoX^lQul"^ /T ^{qi)) versus \Qu\'^, Eq.( |lll|) shows that at small 
q we should find a straight line with intercept 7"^ and slope 7^^A^j^. 

In Figs. ^ — c we show such plots for the three types of perturbations shown in Fig. |T], for 
Nz = 30 and selected values of temperature, fi = y,x, z correspond to the tilt, compression, 
and shear perturbations respectively. The straight lines through the data result from least 
squares fits, using the smallest eight values of g > 0. In virtually all cases, the fit is quite 
reasonable. In Figs. ^ — c and ^ — c we show the values of 7^ and (A^r/A)^ obtained from 
these fits. In each case we show the result of fits to the smallest eight, seven, six, and five 
values of g > 0. As is seen, our results are virtually insensitive to the number of values of q 
used, except for the case of the compression perturbation in the vicinity of Tm ~ 1.2, where 
our data is rather scattered and statistical errors are large (see data for T = 1.2 in Fig. ^6 
and the corresponding dashed line fit). We have also obtained values of 7^ and (A^r/A)^ by 
fitting JoA^|(5;/p/T^(gz/) to a second order polynomial in |Qyp. We have found the results 
from such quadratic fits to be essentially unchanged from the values obtained from the linear 
fits. 

In Figs.|^ — 6 we show 7^ and (A^r/A)^ respectively for all three types of perturbations, 
comparing the finite size effects for Nz = 30 and A^^ = 15. The results shown are for fits to 
the smallest eight values of q, except for the case of the tilt perturbation for A^^ = 15 where 
we have used only the smallest four values of q (since the allowed values of q^ are spaced 
twice as far apart for N^ = 15 as compared to Nz = 30). We see that finite size effects are 
in general small, except for the case of the shear perturbation /i = z. 

We now discuss our results for 7^. From Fig.|^ we see that all three 7^ ~ 1 at low 
temperatures. For 7^, this is in agreement with our expectation that there is a total Meissner 
effect for shear perturbations in the vortex line lattice phase. However the elastic theory 
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results given by the first lines of Eqs.(^) and (|69|) would lead one to expect 7y, |7x| -C 1. 
This is because for the relatively large B simulated here, B c::^ H and the susceptibilities 
dBcr/dHcr that enter 7^ and 'jx in Eqs.(pO|) and (^) are both close to unity. That we 
find 7j;,7x — 1 at low T is, we believe, an artifact of our discretizing grid which acts like 
a periodic pinning potential for vortex lines. At low T, the vortex lines are locked into 
a lattice structure commensurate with this pinning potential. Indeed, the fact that the 
ground state of Fig.|| is not a perfect equilateral triangular lattice is due to this effect. 
This periodic pinning potential leads to an enhanced stiffness of the effective elastic moduli, 
greatly reducing the susceptibilities dBcr/dHcr from their continuum values, and resulting in 
the observed 7j,,7x — 1 at low T. Indeed the periodic pinning potential of the discretizing 
grid acts in many ways like the columnar pins of the "Bose glass" modeled of a disordered 
superconductor, and 7^ = 1 is similar to the "transverse Meissner" effect for tilting the 
applied magnetic field that is found in that problem. One can wonder whether the decrease 
of jy from unity which begins at T ~ 0.6 is a smooth cross-over due to finite energy barriers 
for discretized vortex fiuctuations, or is rather a sharp transition, being the periodic pinning 
analog of the Bose glass transition. 

At higher temperatures, 'j^ and 7^ decrease towards zero at Tm — 1.2. We will soon 
see that this Tm is the melting temperature of the vortex line lattice. It is only when the 
vortex lattice melts that the vortex lines also depin from the the periodic potential of the 
grid. Assuming that the effective tilt and compression moduli of the unpinned vortex line 
liquid are not greatly different from those of the continuum vortex line lattice, one expects 
from Eq.(|T]), < 7j, ^ l/(87rA26o) = 0.024 < 1, and from Eq.(|6|), 7^ ^ -l/(167rA26o) = 
—0.012 < 0. Looking at Figs.|^a — 6 we see that 7^ is indeed small and positive for T > T^, 
while 7^. is small and negative. That 'j^ is indeed negative and not zero for T > T^ is 
more clearly seen in Fig.^ by noting that the intercepts of the fitted lines are just 7"^. 
The numerical values we find for 7^; and jy in this region are in roughly the same ratio 
as the above estimates, but approximately two or three times larger in magnitude. This 
rough agreement of 'jx and 7^ with elastic theory gives us confidence that, above Tm, the 
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artificial pinning introduced by our discretizing grid is no longer playing a significant role in 
the vortex line fiuctuations. 

Returning to Fig.|^a we see that, in contrast to jy and •jx, Iz remains equal to unity well 
into the vortex line liquid phase T > T^. 7^ only decreases from unity towards the small 
value expected from Eq.(P5|) for the vortex line liquid, 7^ = 1 — dB±/dH± = 7^, at T^ ^ 1.8. 
Thus longitudinal superconductivity, indicated by the shear Meissner effect with 7^ = 1, 
persists well above Tm into the vortex line liquid. This one of the main observations of our 
simulations. Comparing results for A^^ = 15 with Nz = 30, we see that the width of this 
transition clearly sharpens as Nz increases, however the temperature Tc, where 72 starts to 
fall below unity, decreases only slightly. 

We now consider our results for (A^r/A)^. For the tilt perturbation, comparison of 
Figs. |a and |]a show that to a very good approximation, 7^ ^ (Ay^j/A)^ for the entire range 
of T. Such a result follows from Eqs.(^) and ( pSD if one makes the simple assumption 
that C44(q) ~ {B'^/4:nX\)Vg± + boii where Vg± is the vortex line interaction of Eq. (pTj) , and 
ei = 77^^ei where ei is the effective Qz independent single vortex line tension. It is interesting 
that Xy^ shows no increase as Tm is approached from below, as is usually associated with a 
decay length near a transition. 

Turning to the compression perturbation we see from Fig. ^ that, in contrast to A^r, 
(Axr/A)^ does increase from unity as T^ is approached from below. This increase is clearly 
noticeable at temperatures sufficiently below T^ so that our data still has good statistical 
accuracy. This is in contrast to a similar increase in 73. in Fig.|^ just below T^, which we 
do not believe is statistically meaningful, but is rather just a refiection of the scatter in our 
data. Thus, as the lattice starts to depin from the discretizing grid, a fiuctuation in vortex 
line density decays over an increasing length scale X^r- Above T^, our numerical values 
are consistent with (A^r/A)^ ~ •jx, as expected from Eqs. (|69D and (^) for the case A > a^. 
That A^R is indeed negative here, and so Xxr is imaginary, can be seen in Fig. ^ by noting 
that the slopes of the fitted lines are equal to 7;^^A^j^, and that for T > T-^, 7^. < 0. 

Finally, we turn to the shear perturbation. Since this perturbation experiences a total 
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Meissner screening in the superconducting state, we may expect, in analogy with the Meiss- 
ner effect at H = 0, that A^j^ ~ Ug where Ug is the density of superconducting electrons 
(not to be confused with ps boson, the superfiuid density of the analog 2D bosons). Since rig 
decreases as T increases, vanishing at the superconducting transition, we expect that A^r 
should increase with increasing T, reaching a maximum at T^ (diverging in the case of a 
second order transition). Precisely such behavior is seen in Fig. ^. Above Tc, A^r decreases 
to roughly the same small values as Aj^r, as is expected from Eqs.(^5|) and (|8^ . Comparing 
results for A^^ = 15 with A^^ = 30 we see that, similar to the behavior of 7^, the transition 
width sharpens and the temperature of the peak in A^r slightly decreases as Nz is increased. 
It is interesting to note however, that the value of A^r at its peak has also very slightly 
decreased as N^ increased. 

The possibility that longitudinal superconductivity can persist into the vortex line liquid 
has been suggested by the 2D boson analogy. We can therefore compare the Tc found here 
with the predictions of Sec. |V B| . Rewriting Eq. (|86|) in terms of the dimensionless parameters 



of our numerical simulation, and taking dHj_/dB_j_ ^ 1, gives T^ = %tiJqk^ /N~,. Using k = 5 
and Nz = 30 gives T^/Jq ~ 21, ten times larger than the value 1.8 found numerically. We 
can also estimate the entanglement cross-over of Nelson. Eq.(P^) gives Tx = ttJo In K/2boNz, 
and using bo = 1/15, k = 5, and Nz = 30 gives Tx = 1.26 ~ T^. This is somewhat lower 
than the observed T^. Moreover, both the boson superfiuid transition temperature and 
the entanglement temperature Tx should scale with system thickness as l/Nz- In contrast, 
comparing A^^ = 15 with A^^ = 30, we see no such dramatic shift in the numerically observed 
Tc ~ 1.8. 

D. Results: Vortex Line Fluctuations 

To elucidate the nature of the transitions in our model, we have measured other properties 
to characterize the vortex line fluctuations in the system. In Fig. ^ we show snapshot views 
of the vortex line configurations for A'^ = 15, at various temperatures T < T^, T,^ < T < T^ 
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and Tc < T. We show both a side perspective and a view looking down along the applied 
field. We see clearly that for T < T^ there is a vortex line lattice. For T^ < T < T^ the 
lattice is disordered but the vortex lines remain for the most part disentangled. For Tc < T 
the lines are highly entangled. 

For a quantitative determination of the vortex line lattice melting temperature, we com- 
pute the structure function of vortices within the same xy plane, 

Si^±) = ^T.e'''^'^'''"''Hn^.n,.)S.^,.^ (116) 

Below Tm we expect to see Bragg peaks at the reciprocal lattice vectors K of the vortex 
line lattice, while above Tm we expect to see approximately circular rings characteristic of 
a liquid. Let us denote by {Ki} the six smallest non-zero reciprocal lattice vectors, and by 
{K'^} the six vectors obtained by reflecting the {Ki} through the x axis. Then since the 
ground state vortex lattice of Fig. ^ breaks this reflection symmetry, while the vortex line 
liquid restores it, the quantity AS = ^(Ki) — ^(K'^), averaged over the six {Ki}, serves 
as a convenient order parameter for the melting transition. We plot AS*, normalized by 
5*0 = ^(K = 0), in Fig.|^. We see that AS* vanishes at T^ — 1.2. In an earlier workM we 
have shown intensity plots of 5'(q_L) in the entire q_L plane. The circular rings seen above 
Tm verify that T^ is indeed a melting to a liquid, and not a depinning to a floating vortex 
lattice, or some other vortex lattice structural transition. 

As another measure of vortex line fluctuations, we have computed the fluctuation length 
of the vortex lines in the directions transverse and parallel to the applied magnetic fleld. 
The total length of vortex lines in the ground state is Cq = bQNzN]_. If, in any conflguration, 
Cf^ is the total length of all vortex lines in direction fi (we count length here as an absolute 
quantity; oppositely oriented segments do not cancel each other out), then we deflne the 
normalized excess vortex line lengths as Ai± = {C^ + ^y)/{'^^o) and Af^ = {Cz — Co)/Co- 



We plot Ai± and Ai^ in Fig.jiy. If we assume that all vortex fluctuations consist of purely 
transverse motion of the magnetic fleld induced lines, then Ai± is the average transverse 
distance traveled by a vortex line between two adjacent xy planes. If we further assume that 
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these lines are fluctuating as in a random walk, then the total transverse deflection of a line 
in traveling down the entire length of the system A^^ is, u = \^N^A£±. Entanglement should 
occur when u ~ a^, or when Ai± ^ a^/^/Wz. From Fig. ^, and using a^ ^ l/-\/5o = 3.87, 
we would estimate the entanglement temperatures as Tx ~ 2.1 for N^ = 15, and Tx ^ 1.9 
for A^^ = 30. These are both consistent with the Tc seen in Fig.|^. However, if the transition 
at Tc is indeed caused by the onset of entanglement due to transverse wandering of magnetic 
fleld induced vortex lines, it is necessary to explain how just above Tc, where 7^ ~ 0, one 
can have a Tboson/Tboson — NzT /{At^'^Jqk^) ^ 0.06 (see Eq. (|43|) ) so much smaller than the 
lower bound 2/7r given by the Kosterlitz-Thouless theoryEj of the analog boson superfluid 
transition. We further note that in previous simulationgHZl'ca with A — ^ cxo, where samples 
up to thickness A^^ = 200 were studied, the above criterion gives a Tx which is well below 
the observed Tp. 



Returning now to Fig. |T0|, we see that the above assumption of strictly transverse fluctua- 
tions of the fleld induced lines, while reasonable near the melting T^ where A^^/Af^ — 0.035, 
is not at all reasonable near Tc-, where A^^/Ai*^ ~ 0.41. The excess vorticity along z can 
only come from either fleld induced lines which wander backwards, or from closed vortex 
ring excitations. Both these types of excitations are absent from the usual 2D boson analogy. 

Using an algorithm we have describe elsewhere,!^ we trace out the vortex line paths in our 
conflgurations to compute the distribution q{p) of the number of closed rings of perimeter 
p, normalized by the ground state vortex line length Cq = h^NzN^. In Fig.|ll] we compare 
the total length of all vortex line fluctuations, Aftot = 2A£^ + Af^, with the total length 
of all vortex ring excitations, A^ing = Y.pVlip)- We see that A^ing ^ A£tot through the 



melting Tm, however at Tc = 1.8, A^ing has increased to 27% of A^tot- In Fig.|T2| we show a 
semi-log plot of q{p) vs. 1/T. The straight lines found at low T indicate thermal activation 
with a constant energy barrier that increases with ring size. At high T ~ 2.8 the q{p) curves 
saturate. Note that the thermal activation for rings persists up to temperatures above Tc. 
This suggests that, although the number of rings is becoming sizable near Tc, the transition 
at Tc is not directly associated with any critical behavior of the rings. This behavior is the 

42 



same that we saw in simulations of a 3D XY model, corresponding to A ^ oo, when we took 
anisotropic couplings;Bj for isotropic couplingsE3 in the XY model, the saturation of the q{p) 
curves coincided with Tc. In Fig. |1^ we plot the specific heat C vs. T, for A^^ = 30. We 
see that C rises smoothly through Tc. The peak occurs near T ~ 3.0 (we only have enough 
data at high T to locate it very crudely), where the q{p) curves saturate. The peak in C is 
thus associated with the proliferation of the closed vortex rings, which we believe to be a 
non-singular cross-over phenomenon associated with the transition of the zero field 60 = 
model, which occursE^I at Tco ~ 3. The peak in C is also probably associatediHa with the 
onset of a strong diamagnetic response in the system, which occurs at the so-called "mean 
field H,2{Ty line. 

Finally, we consider the entanglement of the magnetic field induced vortex lines. Due to 
the periodic boundary conditions along z, the set of points {r_ii{Nz)} where the field induced 
vortex lines pierce the xy plane at ^ = iV^, must be some permutation of the set of points 
{r_Li(0)} where the lines pierce the xy plane at z = 0. Lines for which rj_i{Nz) = r±j{0), 
with i ^ j, form part of an entangled braid when viewed in the periodically repeated system. 
We can thus classify each magnetic field induced line as belonging to a given braid of order 
m, according to the number of lines m that are mutually entangled in the preceding sense. 
We compute the distribution n{m) giving the average number of lines in a braid of order m, 
where Y^m^i^ = ^o^l is just the total number of field induced lines. 

In the 2D boson analogy, such entangled vortex lines represent particle exchanges. A 
superfluid state of these 2D bosons is expected when there are many such exchanges, and in 
particular when there is a finite probability to form large exchanges involving a macroscopic 
fraction of the particles,E3 that wrap entirely around the system in the transverse direction 
and thus contribute to 12^0 = liniq^o('^y(Q'^)^j;(— gx)). 7^, is a direct measure of n^o (see 
Eq.(p8D), and hence a measure of the presence of such large exchanges. In Fig.|T^ we plot 
vs. T the fraction of lines R = n{l)/boN]_ which are not involved in any particle exchanges, 
i.e. the fraction of unentangled vortex lines for which r±i{Nz) = r_Lj(0). We see that -R = 1 
and all lines remain unentangled up to T ^ Tc, at which point R decreases towards zero. 
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The width of the decrease in R is roughly the same as the width of the decrease in 7^, for 
both A^, = 15 and 30. 



In Fig. ^ we plot the entanglement distribution n{m) vs. m, for several values of T near 
and above Tc = 1.8, for A''^ = 30. We see that the distribution broadens as T increases, 
indicating larger particle exchanges, however no sharp feature is obvious as T increases 
through Tc. This is in contrast to what we observed in simulationscZl'Ea of the A — i> cxo 
3D XY model, where n{m) got dramatically flat and equal to unity over a wide range of 
intermediate m as T reached Tc; these XY simulations however used much larger system 
thicknesses, A^^ ~ 100 — 200, and this is one possible reason for the difference in behavior 
from the present case. 

VII. CONCLUSIONS AND DISCUSSION 

The main conclusion of our numerical work is that longitudinal superconductivity van- 
ishes at a Tc which lies well within the vortex line liquid, at least for the system sizes we have 
been able to investigate. We note that our system sizes N± = 30, and A^^ = 15, 30, are large 
compared to the microscopic length scales of our model, A/^_l = 5 and a^/^±_ — vT5. We 
have discussed a mechanism for this phenomenon in terms of the KT superfiuid transition of 
the analog 2D bosons. However the Tc predicted by Eg . (|86|) is an order of magnitude larger 
than the numerically observed value. The entanglement temperature Tx of Nelson is of the 



correct order of magnitude as the observed Tc. Figs. H and |T4| also suggest a connection 
between geometrical entanglement and Tc. However upon comparing N^ = 15 and 30, we 
failed to see any sign of the dramatic size dependence Tx oc l/N^ that is expected from 
Eq. (p4D . A similar size dependence is also expected for the KT prediction of Eq.(^). 

In earlier simulation^Hea of a 3D XY model, corresponding to the A ^ 00 approximation 



of Sec. [Ill D|, we have studied much thicker systems than reported on here, with N^ as large 



as 200. We again found longitudinal superconductivity to vanish at a T^ within the vortex 
line liquid, with virtually no finite size effects in the apparent value of T^ as A^^ was varied. 
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An analysisEZia of geometrical entanglement, as done here in connection with Fig. 0, gives 
a Tx well below the observed T^ for the thicker systems, and the dependence of T^ on 
the system anisotropy was found^ to be T^ oc 1/r/, rather than the Tx oc l/r^^ predicted 
by Eg. (p^) . New simulations" have further shown that there is no apparent change in the 
large N^ limiting value of T^ when the periodic boundary conditions along the direction z 
of the applied magnetic field are replaced with the more realistic free boundary conditions. 
We believe that these A — i> oo simulations are therefore in good agreement with the work 
of Feigelman and co-workers,l2ilc3 who argued for just such a superconducting to normal 
vortex line liquid transition, with a T^ which remains finite as L^ — > cxd. We believe that 
this transition of Feigelman et al. applies strictly to the A — > oo model, and represents a 
^boson = Meissner transition, as ^boson varies, for the analog 2D charged bosons. 

Returning to our present simulations, we believe that our results represent a finite A 
cross-over from the above A ^ oo transition at T^ . Although we believe that the A — i> oo 
limit is extremely subtle, one may imagine the following scenario. When A is large, although 
the analog 2D bosons behave like a neutral superfiuid on sufficiently long transverse length 
scales, on small length scales they will have the A — i> oo behavior of charged bosons. We 
would then expect the 2D boson hehcity modulus to have, at finite transverse wavevector g, 
a piece that looks like that of Eq.(|50|). We thus expect a form like, 

[<^A ] boson Q' 

^1 _|_ \2 ~2 

^ "T ^Rbosony 



-*- boson l^yj J-bosonl^UJ ~r /boson -, i\2 2' 1/ 



As g -^ 0, it is Tboson(O) that determines if the 2D bosons are in a superfiuid (Tboson(O) > 0) 
or a normal fiuid (Tboson(O) = 0) state; but at sufficiently large q it will be the second term 
that dominates, giving the appearance of a charged boson system. For T < Tc, with Tc 
the 2D neutral boson superfiuid transition of Eq. (p6|) , one has Tboson(O) = and only the 
second term is present. As g — > 0, this term vanishes, and Eq.p3|) then gives 7^ = 1, i.e. 
we have the perfect shear Meissner effect that we expect for the 2D boson normal fluid 
phase, as discussed in Sec. [Ill C|. However, if L^ is thin enough that T^ <^ T^ with T^ the 



Meissner transition of the A ^ 00 2D charged boson model, then as one cools down to T^, 
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one expects ARboson will become large, and possibly of order the finite transverse size L^ of 
the system. In this case, for all available wavevectors, q > 2tx /L±_ yields y^\\,osoa'f' ^ 1' ^^^ 
the second term becomes approximately the constant, 7boson[-'^A^]boson/ARboson- Eq.(^3]) then 
gives 72 = 1 — [7bosonA^/A|^tiogo^]. It thus appears as if the perfect shear Meissner effect has 
been lost at the lower temperature ~ T^ . We note that for this scenario to agree with the 
small values of 72 that we find at temperatures above our numerically observed value of Tc, 
it would be necessary to have 7bosonA5_/ARboson ~ Ij i^ is not apriori obvious why this would 
be so. 

Thus for a finite A simulation to see other than the above A ^ 00 cross-over behavior, 
it would be necessary to do one of the following. One could increase the transverse size L^, 
keeping L^ constant, until one is in the hmit where max[ARboson] ^ L^ (although in the 
A ^ 00 model ARboson might diverge at T^°°, in a finite A model any such divergence would 
be rounded out to a finite maximum value). In this limit, the second term in Eq.( [117| ) would 
be observed to be ~ g^, and so one would find 7^ = 1. One thus expects the apparent T^ to 
increase above T^ as L_|_ increases above max[ARboson]- One might never actually reach the 
true 2D boson neutral superfluid transition T^ of Eq.(|^), since as temperature increases, 
thermally excited closed vortex rings will start to proliferate, and vortex lines can make long 
transverse wanderings between two adjacent xy planes; both such effects are left out of the 
naive mapping to 2D boson statistical mechanics. Alternatively, one could keep Lx_ constant, 
but increase L^, so that the T^ of Eq.(^) falls below T^ . For our parameters, Eq.(^3D 
suggests that this would require a system of thickness A^^ = (J)1/2tt'^TI^ = Sir Jqk'^/T^ ^ 320. 

Several other groups have done simulations similar to ours. Most of theseEHH have been 
in the A — i> 00 limit, but at much higher vortex line densities such as 60 = 1/6. In these 
cases it was found that Tc ~ T^, and so no longitudinal superconductivity was observed in 
the vortex line liquid. We believe that this is a consequence of the high densities Bq which 
have been used. Recently, we have studied" the phase diagram in such A ^ 00 XY models, 
as a function of the system anisotropy rj. Increasing 77 at fixed 60 can be argued to play 
a role similar to increasing bo at fixed rj. We found that as 77 increased, Tc and Tm came 
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closer together, and eventually became indistinguishable from each other. Similar results 
have recently been reported in simulations by Koshelev.0 

Sasik and Stroud have done simulations^ for the A ^ cxd limit using the lowest Landau 
level approximation, which treats the xy planes as a continuum and so avoids the artificial 
pinning of our discretized London model. For all values of anisotropy studied they find 
Tc ~ T^. However TesanovicE^ has argued that the lowest Landau level approximation fails 
as the magnetic field decreases, and so at such low magnetic fields, the London and the 
lowest Landau level approaches need not be in agreement. Using a mean field analysis, 
TesanovicE3 has argued that longitudinal superconductivity can persist into the vortex line 
liquid in this low field limit. 

Finite A simulations have been carried out, for the same discretized London model as 
considered here, by Carneiro.cl'o For large line densities, he finds Tc ~ Tm, consistent with the 
above A — i> oo results. For line densities comparable to our own, he finds Tc noticeably above 
Tm, when following our analysis based on the q dependence of the helicity modulus within 
the Helmholtz ensemble of fixed internal magnetic field 60 • He has suggestedEJ however 
that the result may be very sensitive to the q -^ extrapolation implied by fitting to the 
expansion of Eq.( |110D , with different results obtained when truncating at different orders of 



the expansion, or when using a different number of q data points in the fit. However our 
fits of Figs.^ and ^ show essentially no sensitivity to the number of q data points used, or 
when comparing a linear versus a quadratic order fit to the data of Fig. ^. Carneiro has also 
carried out simulationsoo in a Gibbs ensemble, in which the total transverse magnetic field is 
allowed to fluctuate. Here he concludes Tc ~ Tm, even for dilute densities comparable to our 
own. However we believer that in this case, his q = calculation of the fluctuation in the 
transverse magnetic field cannot distinguish between the shear perturbation of Fig.|I|c, which 
is related to the 2D analog boson superfluid density, and the tilt perturbation of Fig. |l]a, 
which is not. We believe that his results are reflecting the softening of C44 that occurs at the 
depinning/melting transition (as is observed in our Fig.|^), rather than reflecting the loss 
of longitudinal superconductivity. 
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Most recently, simulations at finite A have been carried out by Nguyen et a/.,E^ who 
extend our work to consider behavior as the anisotropy rj is varied. For an isotropic system, 
they find Tc well above Tm, in good agreement with our results. However as rj increases, 
they find the very intriguing result that Tc decreases, and eventually falls below T^. Such a 
possibility (not observed in similar A — i> cxo simulationg^) has been suggested by Frey et alB^ 
as a result of dislocations proliferating in the vortex line lattice. Glazman and KoshelevH 
have made similar predictions, based on the effect that vortex lattice elastic fluctuations 
have in reducing the interplanar Josephson coupling. However Nguyen et al. suggest that 
their result is due to the proliferation of vortex rings between adjacent xy plans, and they 
flnd at high anisotropy that T^ ~ r^"^, rather than the T^ ~ rj^^ predicted by Ref. [ ^ or 
the Tc ~ l/ln?7 predicted by Ref.[ 0. It should be noted however that Nguyen el al. 
base their criterion for superconductivity on computing the helicity modulus at the single 
smallest non-zero value of q allowed by their flnite size system. We have argued above that 
a more careful analysis should be based on parameters extracted from the q dependence of 
the helicity modulus, as g ^ 0. Conclusions based on calculations at speciflc values of flnite 
q can more easily be led astray by subtle cross-over effects such as we have discussed above. 
Clearly more systematic studies, using our g ^ analysis, and making a more extensive 
study of flnite size dependencies, need to be done for both the isotropic and anisotropic 
cases. 

Our result of Eq.(p9D suggests that one should flnd T^ > T^, and hence longitudinal 
superconductivity within a region of the vortex line liquid, whenever a sample is thinner 
than Lzraax ~ 1400/im, for a Tm ~ 90°K. Virtually all experimental single crystal samples 
fall below this critical thickness. One can therefore ask whether any experimental evidence 
favors our conclusions. Naively, one would expect a vortex line liquid with longitudinal 
superconductivity to show a flnite linear resistivity transverse to the applied magnetic fleld, 
but zero linear resistivity parallel to the applied fleld. However, in MC simulations of a 
A ^ oo XY model,c3 we found that in the intermediate phase Tm < T < Tc vortex density 
correlations decayed anomalously slowly with time. This suggests that vortex lines may 
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be moving more slowly than diffusion, and if so, it is not obvious what to expect for the 
transverse resistivity. Experimentally, it is very difficult to obtain accurate measurements 
of the longitudinal resistivity, due to the slab geometry of single crystal samples, and the 
non- uniformity of current distributions. Transverse resistivity measurements are intimately 
related to the vortex pinning impurities in the sample, and so are also not unambiguously 
characterized. 

Nevertheless, the following suggestive observations have been made. Experiments by 
Steel et alM on artificially prepared MoGe/Ge layered superconductors, found that the d.c. 
resistivity parallel to the applied field decreased sharply, and showed an onset of strong 
nonlinear behavior, at a temperature above that where the transverse resistivity vanished. 
Experiments by Kwok et alE3 on YBCO, studying the pinning of vortex lines to twin grain 
boundaries in a system with a well controlled small number of twin planes, found evidence 
for a sharp lock-in pinning transition at a temperature above vortex lattice melting (where 
the melting transition was determined by the observation of a sharp drop in transverse 
resistivity). Such a lock-in transition within the vortex line liquid may suggest a transition in 
the nature of vortex line fluctuations, as at our T^. Early experiments by Safar et al.f^ using 
a flux transformer geometry, similarly showed evidence for the onset of coherence parallel to 
the applied field at a "Tt/i" above the temperature where the transverse resistivity vanished. 
However more recent flux transformer experiments by Lopez et a/.H showed that these two 
temperatures merged when the sample was made purer, with all twin grain boundaries 
eliminated. Moorecll has recently proposed an interpretation which argues that the single 
transition observed in these newer transformer experiments is the result of some very rapidly 
increasing longitudinal length scale, rather than being a first order vortex lattice melting 
transition, as is the usual interpretation. If correct, such a rapidly increasing longitudinal 
length scale might be associated with our T^. 

The experimental evidence cited above remains, at best, inconclusive. There are several 
possible reasons why observing a T^ > Tm might be experimentally difficult. Firstly, as 
in our simulations, the relevant temperature is likely to be the A — i> oo transition T^, 
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rather than the much higher T^ of Eg. (|86|) . In recent simulationsE^ of the A ^ oo XY 
model we found that T^ and Tm merged as the anisotropy 77 increased. How far apart the 
corresponding T^ and T^ for any particular real material are likely to be, remains unknown. 
Secondly, real layered high temperature superconductors are likely to have an interplanar 
Josephson coupling that is proportional to the cosine of the phase angle difference across 
adjacent planes. The non quadratic nature of such a cosine interaction leads to a coupling 
between spin wave and vortex fluctuations that is absent in both our continuum model and 
our discretized model of Eq. (|97|) using the Villain interaction. As either magnetic field, 
temperature, or anisotropy increases, large interplanar phase differences can be induced by 
elastic vortex line fluctuations, leading to a large decrease in the effective interplanar energy 
coupling constant. Such a "decoupling" cross-over, as discussed by Glazman and Koshelev,El 
and Daemen et al.fn might obscure any true critical behavior at a higher Tc- Finally, the 
free boundary conditions of a real superconductor, as opposed to the periodic boundary 
conditions of the 2D boson mapping and of our simulations, might lead to a more effective 
washing out of the 2D boson superfiuid transitioncHl than we have imagined, for samples of 
experimental thickness. 
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APPENDIX A: SUPERFLUID DENSITY OF 2D BOSONS IN THE PATH 

INTEGRAL FORMULATION 

The superfiuid density of a system of A^ interacting bosons can be defined in terms of the 
response of the system to the presence of a heat bath moving with velocity v(r) ("moving 
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walls"). In the following, all position, velocity, and wave vectors are two dimensional vectors 
in the xy plane. 

The average momentum density {pq^)v that results in linear response to the heat bath 
velocity V-qy = jj^ J d'^rv^{r)e~'''^'^ is given by 

{Pqt^)v = x,^u{q)v-q^. (Al) 

For an isotropic 2D system, the momentum density susceptibility can be written in terms 
of its longitudinal and transverse pieces, 

Xfiuiq) = qA'^XLiq) + [Sf^u - qMxAq)- (A2) 

The number density of superfluid bosons ps is then given in terms of the transverse suscep- 
tibility by§'§ 

mps = mp — \iTCiXTiq) (A3) 

where m is the boson particle mass, and p = N/L"^ is the total boson density. 

For a system of interacting bosons in the presence of a moving heat bath, the Hamiltonian 
in the reference frame of the heat bath is given by, 

^ = E ^iP^ - "^v(r,))2 + V({r, - r,}) (A4) 

i 

where the interaction V depends only on the bosons relative positions. 

The partition function is 2 = Tr[e^'^^], and the free energy is JF = —Tin 2, with 
T = 1//5. Consider now that v points only in the y direction and varies only in the x 
direction (so the v is purely transverse). If we write ?i[v] = ?-^[0] + 5?i[v], then since in the 
g — i> limit that v becomes a uniform constant 67i and 7Y[0] commute, one has 

q^o dvy{qit)dvy{—qx) 
To evaluate ps in terms of the path integral formalism,!^ one now writes the Lagrangian 
associated with the Hamiltonian of Eg. ( |A4D , and transforms from real time t to imaginary 
time r = it. One gets 
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mp — limxT(gx) = mps (A5) 



ti=0 



A^) = -Ey(^) -n{r.-r,}) + zmE(^^).v(r.). (A6) 

The partition function is then given by 

Z = /'l?[{ri(r)}]e^"'/o"''^-^W. (A7) 

where the sum is over all possible boson world lines {rj(r)} subject to permuted periodic 
boundary conditions, i.e. {rj(0)} = 7'{rj(/5)} where V is any permutation of the A^ bosons. 
Applying Eq. ([A5| ) to the above form for Z then results in 



mps = lim 



9-0 L'^n 



2 



(r''^i:^^-)(r*'i:^--))^ m 



where (. . .)o denotes an average over world lines weighted by the Lagrangian factor as in 
Eg . ([A^) , only now taking v = in £. 



Note that the heat bath velocity v(r) enters the Hamiltonian (|A4|) and the Lagrangian 
with precisely the same form as would a 2D external magnetic vector potential given by 
V = {h/m)A'^^^ (where, as in Sec.|I|, the units of A'^''* are such that V x A*'''* = (27r/0o)-f^z, 
with ffz the 2D magnetic field). In analogy with Eq.(^) we can thus define the helicity 
modulus of the 2D bosons as, 

i boson (Q') = -^ r^ 4 / — ^TTTT"? ^ ~ 



dAy{qx.)dAy{—qx) m? dvy{qx)dvy{—qx.) 

1? 



r^rZ'-^e-«)(f,r'Z'i^e'-'A) , (A9) 



with. 



limTboson(g) = -Ps = T{W^)o , (AlO) 

where Wy is the y component of the "winding number" introduced by Pollock and CeperleyES 
in their path integral approach to the superfluid transition in boson systems. 

Our derivation above can be modified in a straightforward way to deal with a boson 
interaction mediated by a gauge field, as is the case for the more realistic London interaction 
between vortex lines.E^ One just replaces the pair potential V^({ri — Tj}) with the necessary 
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coupling to the gauge field, and free field energy terms. However the coupling of the bosons 
to an external vector potential remains unchanged. Thus the expression for the 2D boson 
helicity modulus in terms of boson world lines remains unchanged from Eq. 

APPENDIX B: ELASTIC MODULI 

In this appendix we summarize some results concerning the elastic moduli which ap- 
pear in Eq . (|5^ . Although calculations of these moduli have appeared elsewhere Jl^'EJo our 
explicit computation of the order q^ dependence at small q we believe is new. 

As shown by Sudb0 and Brandt ,c3 the elastic tensor $Q,^(q) can be expressed in terms 
of the vortex line interaction tensor Kt/3(q), as 

^aM = -^ E {?'^«/3(K - q) + (K - q)„(K - Ci)pV,,{V. - q) - K^KpV,,{K)} (Bl) 

where {K} are the reciprocal lattice vectors of the vortex lattice. 

For B = Bti the elastic moduli we are interested in can be expressed in terms of $Q,^(q) 
as 

cm{qj) = —^xx{qy), cii(gx) = — <|)^^(gx), C44(gz) = — $^^(gz). (B2) 



For the London interaction, the sum over K in Eq. (|Bl|) is divergent, and some method 
must be employed to make it converge. As shown by Brandt ,t3 this can be achieved for 
Cm{,<iy) and cii(gx) by subtracting off the self energy of a line interacting with itself. This 
then gives, 

-a.toy) = j^{i:f«lK.«|-f/jg5Fa.|k.,]} (B3) 

Cnfe*)^^{Efn|K.,l-|/|iF„|k.,]}. (B4) 

where we find, after expanding (^^x to O(g^), averaging over the orientation of the vortex 
lattice in the xy plane, and substituting in for Vzz from Eq.(|T^) 
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+ 



(B5) 



+ 



;i + AiF)3 (l + AiA;2)4 (1 + AiP) 



AigM (B6) 



and 



d 



Fn[k, g] = ;Tp {A;V,,(k) + |A;V,,(k) + [fcV,,(k) + |fcV,,(k) + ^fc'V^..(k)] g^} (B7) 



^A^ 



2 ^ [ (1 + AiA;2)2 (l + AiA;2)3 

1 
+ 



9 



+ 



10 



^W) (B8) 



where Vzz = dV^z/dk'^. 

To treat the tih modulus C44(gz), self interactions of the vortex lines are important. One 
therefore handles the convergence of the sum in Eq.( |BlD by introducing a convergence factor 
into the London interaction of Eg. (|BTD , Vgais — > V^^^ = g{C,±q± + il(fz)yqap- Here g{x) — »• 1 
for X < 1, (^(x) — > for x > 1, and one uses an anisotropic cutoff to model the vortex core, 
C,z/^± = A^/A^ = l/rj. Averaging over the orientation of the vortex lattice in the xy plane, 
and substituting in for V"^^ from Eg . ([161) we find 



cu[qz) 



AttX] 



E^44[K,g] 



(B9) 



K 



where 



F44[k, q] = {v:M + ifc'K';(k) + [^^.(k) + \k'Kli^)] q'} 



xlg 



+ 



xlg 



2 [1 + A2P ' (l + AiA;2)2 
Ai^ , Ai^ 



+ g 



g 



+ 



1 + AiP 
g g 



-? + 



g 



1 + A2p)2 (l + AiP)3 1 + A2A;2 (l + AiA;2)2 2 2(1 + A^F 



where V/;^ = dV^^^/dq^ and g = dg/dq^ 



ll^l 



We consider first the limit of large magnetic fields, A_l ^ a^ {a^ is the spacing between 
vortex lines). In this case one can approximate the sum over K by 



(BIO) 



Alg2 (Bir 
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J2F[K,q] = F[0,q] + ^^^^J^^ dkkF[k,q], 



(B12) 



where {AK)"^ = 47r^i?/0o = ^^o ^^ ^^^ area per reciprocal lattice vector, and kg ~ 1/a^ is 
the edge of an approximate circular Brillouin Zone. Carrying out the integrations, we get 



C66(gy) 
cii(gx) 



A_|_/Cq 



52 r. 1 

47r 



A_|_/Co 



1- 



4(1 + Xlki) 



4(1 + Aifcg)4 
3 

" 4(1 + xikiy 



xW 



4(1 + Xlk^f 4(1 + Xlki)^ 



2 „2 



Alg 



j^2 ( 

C44(gz) = — 1 



2AiA;2 



2^2' 



r] ^In 



1 



1 + K^'r/ 



+ 



A_|_K,o 



7] 



-2 



+ 



2 „2 



Aig 



2A^A:g(l + A2fcg) 4Aifcg(l + Aifco2)2 
where for C44 we have taken the cutoff ^_l -^ in all non-divergent terms, and k = A^/^ 
Expanding for large A^fco; we get to the lowest non trivial order 

2 XW \ 



C66(gy) = -r 



cii(gx) 



C44(gzj 



47r 4AiA;g 

47r 

47r 



1 



A2 i-2 



AlA;o^ 



4AiA;2 



1 + 



4AiA;6 



2Aifc2 



V 



In 



1 + 



i/, 



c2 



5 

1 



Aig^ 



1 



(B16) 
(B17) 



Xlkl 



+ 



2 J2 



Xiq 



(B18) 



2A4fc4 4Alfc6_ 

where Hc2 = 0o/4vr^^. 

Next we consider the case of small magnetic fields, A^ <^ a^. Here it is convenient to use 

cPk 



K ^ R 



where 



F[r,g] 



-jkr 



(2vr)^ 



F[k,g] 



(B19) 



and {R} are the direct Bravais lattice vectors of the vortex lattice. 

For the shear and compression moduli, Cge and Cn, the subtraction terms in Eqs. ( |B3D 
and ( P^ ) cause the R = term of the sum in Eq. ( [B19| ) to vanish. Since the range of the 
interaction Vzz is A_l ^ |R|, it will be a good approximation in the sum over R to keep only 
the six smallest vectors with |R| = a^. The Fourier transforms of Eqs. (|B(j|) and (|B(j|), can 
now be obtained with the help of 



(B13) 
(B14) 



(B15) 
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where K,y is the modified Bessel function of the second kind of order z/, whose asymptotic 
form at large x is 

^ (B21) 



K,y{x) ~ W— e 



Keeping only the leading terms in a^/X±, we find 
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64v/2^Ai 

For the tilt modulus C44 there are two cases to consider, depending on the strength of 
the anisotropy. For very small magnetic fields such that A^ ^ a^, all terms in Eq.( [BllD 
may be treated according to the approximation implied by Eq.( |B19|) . Here the R = term 
dominates all others, and we find 
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C44(gz) = T-^TTp {^'^ [2 In (^z*) - 1] + 1 
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>^W} 



(B24) 



For the intermediate case X± <^ a^ <^ X^, we must combine approximations, using Eg. ( pi 2]) 
for terms involving X^k"^, and Eq. (pi9| ) for terms involving A^A;^. We find 
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FIGURES 
FIG. 1. Schematic representation of three possible perturbations of the external magnetic field: 

(a) the tilt perturbation, (b) the compression perturbation, and (c) the shear perturbation. 

FIG. 2. Labeling conventions for the lattice superconductor. Aif^ are directed outwards from 
site i on the bonds of the direct lattice, bi^ are directed inwards towards the dual site i on the 
bonds of the dual lattice, piercing the plaquettes of the direct lattice as shown. 

FIG. 3. Ground state for vortex line density 6o = 1/15 on a cubic grid. Solid circles indicate 
the locations of the straight vortex lines as they pierce the xy plane. 

FIG. 4. Helicity modulus plotted as JoX^\Qu\'^ /T ^{qi)) vs. \Qu\'^ for various values of T. The 
straight lines are fits to Eq.(lll), and determine the parameters 7^ and (A^r/A)^ for (a) the tilt. 



(6) the compression, and (c) the shear perturbations of Fig. ||. 

FIG. 5. Plots of 7^ vs. T as obtained from the straight line fits of Fig.|^, fitting to the 8, 7, 6, 
and 5 smallest values of q, for (a) the tilt, (b) the compression, and (c) the shear perturbations of 
Fig.[l|. Little sensitivity is seen to the number of values of q used in the fit. 

FIG. 6. Plots of (A^pt/A)^ vs. T as obtained from the straight line fits of Fig.|^, fitting to the 8, 
7, 6, and 5 smallest values of q, for (a) the tilt, (b) the compression, and (c) the shear perturbations 
of Fig. ||. Little sensitivity is seen to the number of values of q used in the fit. 

FIG. 7. Finite size comparison of the parameters (a) 7^ and (6) (A^r/A)^, for the tilt (A), the 
compression (O); ^^^ the shear (<(>) perturbation. Open symbols are data for N^ = 15, while solid 
symbols are data for Nz = 30. 

FIG. 8. Snapshots of vortex line configurations for N^ = 15, for (a) T = 1.0 < Tm, (6) 
Tm < T = 1.6 < Tc, and (c) Tc < T = 2.2. The bottom row is the view looking down along the 
applied magnetic field. 
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FIG. 9. Plot of structure function peak heights, AS(Ki)/So, vs. T for N^ = 15 and 30. 

FIG. 10. Average normalized fluctuation length of vortex lines Aiz and Aij_, parallel and 
transverse to the applied magnetic field, vs. T. We see that Aiz ^ A£± for T < T^, indicating 
that there are only transverse fluctuations of the magnetic fleld induced vortex lines. This is no 
longer true near T^.. Open symbols are for N^ = 15, solid symbols are for N^ = 30, and the solid 
lines are guides to the eye. 

FIG. 11. Average normalized length of all vortex line fluctuations A£tot) and average length in 
closed vortex ring excitations A^nng) vs. T. Open symbols are for Nz = 15, and solid symbols are 
for Nz = 30. 

FIG. 12. Semi-log plot of q{p), the distribution of closed vortex rings of perimeter p, vs. 1/T, 
for several values of p. Straight lines at low T indicate thermal activation. q{p) saturates at a 
temperature above Tc. Open symbols are for A^^^ = 15, and solid symbols are for Nz = 30. 

FIG. 13. Speciflc heat C vs. T for Nz = 30. The peak in C occurs above Tc. 

FIG. 14. Fraction of unentangled magnetic fleld induced vortex lines R vs. T, for Nz = 15 and 
30. Lines start to entangle at Tc. 

FIG. 15. Distribution of entanglement of the magnetic field induced vortex lines. n{m) is the 
number of lines that participate in entanglement braids of order m. Data is shown for several 
temperatures near Tc = 1.8, for A^^ = 30. 
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